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16016 EVALUATION SAND WAVES ESTUARY 


KEY WORDS: Dissection; Ecosystems; Energy; Environmental impact 
statements; Equilibrium; Estuaries; Sand waves; Sediment; Spectral analysis; 
Turnovers 


ABSTRACT: The Bella-Williamson rate sediment turnover-organic content 
sediment (RST-OCS) dissection plane method integrates relevant physical, chemical, 
and biological analyses order estimate the chronic impacts dredging 
estuarine ecosystems. spectral analysis finite Fourier transform estuarine 
sand waves presented which only requires spatial sand wave record that 
relatively easy measure contrast stochastic analyses which require both spatial 
and temporal observations sand wave profiles. The total energy content (or, 
equivalently, the statistical variance the spatial sand wave profile) and the 
equilibrium subrange the wave number spectrum provide readily observable 
measures the RST estuary; economic savings, promise significant. 
addition, uniform power law” for the wave number equilibrium subrange 
presented which offers opportunity obtain sand wave frequency spectrum 
without temporal sand wave profile record through dispersion transformation 
proportional the bed shear velocity. 


REFERENCE: Baliga, R., and Hudspeth, Robert T., “Evaluation Sand Waves 
Estuary,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY2, Proc. 
Paper 16016, February, 1981, pp. 161-178 


16020 TOTAL LOAD TRANSPORT ALLUVIAL CHANNELS 


KEY WORDS: Alluvial channels; Bed load; Canals; Hydraulics; Rivers; 
Sediments; Sediment transport; Shear stress 


ABSTRACT: relationship for computation total load transport alluvial channels 
has been developed analysis extensive data from rivers and canals and from 
laboratory studies. The relationship extremely simple and based 
logicalpremises. The accuracy computation total load from this relationship 
quite good comparison that obtained from some the existing relationships. 


REFERENCE: Raju, Kittur Ranga, Garde, Ramachandra J., and Bhardwaj, Ramesh 
C., “Total Load Transport Alluvial Channels,” Journal the Hydraulics Division, 
ASCE, Vol. 107, No. HY2, Proc. Paper 16020, February, 1981, pp. 179-191 


16012 HAZARD ASSESSMENT ARCH DAM SEISMIC ZONE 


KEY WORDS: Dams (arch); Dynamic loads; Earthquake engineering; 
Earthquakes; Fault; Geology; Seismic properties; Stresses; Structural 
engineering; Vibration 


ABSTRACT: Principal faults the Green Lake Arch Dam region were examined and 
assessment was made potentially damaging earthquakes. earthquake which 
can cause the most severe ground motion the site (Maximum Credible Earthquake) 
and another one which has recurrence interval 100 yr, were selected for use 
design. dynamic analysis using the time-history method was performed, showing that 
vertical contraction joints would partially open for short instants under dynamic 
loading. was concluded that the dam can resist the Maximum Credible Earthquake 
without uncontrolled release water from the reservoir, and that the dam can 
withstand the 100-yr earthquake with little damage. 


REFERENCE: Bowes, Donald E., Sison, Antonio X., and Alan L., “Hazard 
Assessment Arch Dam Seismic Zone,” Journal the Hydraulics Division, 
ASCE, Vol. 107, No. HY2, Proc. Paper 16012, February, 1981, pp. 193-208 


16048 MODELING TRANSVERSE MIXING STREAMS 


KEY WORDS: Coordinates; Discharge; Dispersion; Dyes; 
Mixing; Natural streams; Pollutants; Solutes 


ABSTRACT: model transverse mixing natural streams has been tested using 
dye data from field experiment and has been proven successful. The model uses 
natural coordinate system which follows the curvature the stream and based 
the cumulative discharge concept. the numerical solutions, different assumptions are 
used for the factor diffusion, The best accuracy found when local 
values are used together with dispersion coefficient which can vary 
the downstream direction. The changes the size the dye plume can predicted 
fairly well, even with analytical solution. analysis available field values for 
suggests that the dimensionless dispersion coefficient can represented 
function the width-depth ratio and variable which characterizes the curvature 
the stream. The sinuosity possible variable use. Most the data are distributed 
fairly well along one curve which can used obtain estimates for streams 
with small curvatures. 


REFERENCE: Lau, Lam, and Krishnappan, G., “Modeling Transverse Mixing 
Natural Streams,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY2, 
Proc. Paper 16048, February, 1981, pp. 209-226 


16060 OVERLAND FLOW FROM TIME-DISTRIBUTED RAINFALL 


KEY WORDS: Kinematics; Overland flow; Peak discharge; Rainfall; 
Rainfall-runoff relationships; Thunderstorms; Time factors; Time 
concentration 


ABSTRACT: The influence time-varying rainfall overland flow investigated. 
The kinematic wave equations for turbulent flow (across plane, impermeable surface) 
are solved using time-varying rainfall appropriate for thunderstorms. The peak 
discharge shown function surface length, total precipitation, storm 
duration, and time equilibrium for rainfall constant intensity. For rainfall 


durations equal less than the time-of-concentration, there little difference between 
peak discharges estimated using time-varying rainfall and rainfall constant intensity. 
For rainfall long relative duration, the thunderstorm distribution gives much higher 
peak discharges. 


REFERENCE: Allen T., “Overland Flow from Time-Distributed 
Journal the Hydraulics Division, ASCE, Vol. 107, No. HY2, Proc. Paper 
16060, February, 1981, pp. 227-238 
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EVALUATION SAND WAVES ESTUARY 


INTRODUCTION 


The National Environmental Policy Act (NEPA) establishes continuing en- 
vironmental responsibilities for agencies the Federal government [20, Sec 
One these responsibilities the preparation environmental impact 
statements (EIS) for dredging operations which require Federal action [20, Sec 
comprehensive interdisciplinary methodology for identifying the 
chronic environmental impacts dredging operations estuaries has been 
developed Bella and Williamson (2). methodology for incorporating 
evaluation these environmental impacts into environmental strategy has 
been outlined Bella and Overton (1). The power this diagnostic approach 
toward the preparation and evaluation EIS for estuarine dredging lies 
its relatively simple methodology for integrating disciplinary data which are 
easily measured, into two-dimensional dissection plane from which the chronic 
impacts dredging operations may identified. The Bella-Williamson dissection 
plane approach should result significant economic savings due the elimination 
the collection and analyses irrelevant environmental data from estuaries 
which dredging operations are proposed. Included the Bella-Williamson 
dissection plane methodology are estuaries which sand waves are present 
the benthic interface; and the collection and analyses sand wave data 
spectral methods are presented the context disciplinary element 
the holistic dissection plane approach for diagnosing the environmental impacts 
dredging estuaries. 

The two orthogonal dimensions identified Bella and Williamson (2) for 
their dissection plane are: (1) The rate sediment turnover (RST abscissa); 

Research Asst., Dept. Civ. Engrg. and Ocean Engrg. Programs, Oregon 
State Univ., Corvallis, Oreg. 97331. 

Assoc. Prof., Dept. Civ. Engrg. and Ocean Engrg. Programs, Oregon State Univ., 
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and (2) the organic content the sediment (OCS ordinate). Contours measures 
the various physical, chemical, and biological environmental properties 
estuarine ecosystem may displayed this dissection plane with the RST 
and the OCS the axes. From these contours this dissection plane, the 
environmental impacts from dredging the estuarine ecosystem may estimated 
from known predicted alterations the RST and the OCS. Changes the 
RST may diagnosed from the direct influence the hydrodynamic forces 
sand waves the benthic interface estuary. Accordingly, sand waves 
subjected hydrodynamic forces due oscillatory flows estuary were 
observed order estimate the RST pristine estuary. These observations 


from pristine estuary are also compared with sand waves adjacent developed 
estuary. 


Slough 


FIG. Sand Wave Sampling Stations, Coos Bay Estuary, Oregon 


The hydrodynamic forces which directly influence the RST erodible 
sediment interface result alternating sequence erosion and deposition. 
those instances when this alternating sequence erosion and deposition 
result wave-like deformation the sediment boundary, wave propagation 
theories may employed compute the wave celerity, wave length, 
and wave period, these waves. deterministic wave propagation problems, 
the reciprocal the interval periodicity yields the frequency, which 
the process periodically repeated and provides estimate the RST. Since 
the benthic interfaces most estuaries are definitely not purely sinusoidal 
Strictly periodic, the superposition many waves may employed Fourier 
analyze the instantaneous profile this random interface. addition, the 
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hydraulic processes which result erosion and deposition are not, general, 
deterministic and the methods Fourier analyses may further extended 
describe stochastic process spectral methods. 

order demonstrate application the relatively simple spectral analysis 
for estimating the RST estuary from the migration sand waves, sand 
wave measurements the South Slough estuary, Coos Bay, Oregon were 
analyzed. This estuary located the southern Oregon coast approx 341 


south the Columbia river entrance and 702 north the San Francisco 
Bay (see Fig. 1). 


The evolution the various complex analytical models which describe the 
response deformable benthic interface hydrodynamic forces eloquently 
chronicled Graf (7), Raudkivi (17), and Kennedy (9,10). This evolution follows 
essentially the same three conceptual areas described Kinsman (11) for random 
surface gravity waves, i.e., hydrodynamics, stochastic processes and probability, 
and Fourier analysis. The evolution the hydrodynamic concept began with 
Exner (compare Ref. Ref. 17) who assumed that the bottom was initially 
sinusoid and that the free surface the water was sensibly horizontal. Exner’s 
model capable (at least theoretically) including the effects friction 
channel variable width. The principal contribution this model that 
the instantaneous profile the bed form results faster speed for the crest 
the sinusoid compared the trough and provides means for estimating 
approximate time for this sinusoid become unstable and collapse. 
Experimental verification this model lacking, well physical explanation 
how the bed initially becomes sinusoid what shape assumes following 
the instability breaking. 

Others (6,12,15,18,19,21) have extended this model include steady flow 
over sinusoidal bottom with corresponding sinusoidal free surface. The 
linearized potential theory solution derived these authors yields disturbing 
singularity for the special case which the speed the free surface sinusoidal 
profile exactly equal the free stream velocity. The existence this singularity 
frequently used argue that deformable bedform may not remain flat. 
Kennedy (9,10) and Reynolds (18) also have included slowly moving sinusoidal 
bottom their model. However, the speed this motion considered 
small compared with the speed the free-stream velocity and neglected 
the mathematical solution. 

Mei (14) analyzed steady flow with periodic free surface over impermeable 
stationary bottom which also has periodic profile. considering the fully 
nonlinear boundary conditions and assuming that the amplitude the periodic 
bottom profile was higher order the perturbation parameter, successfully 
removed the singularity present the previous linearized studies and was able 
determine the conditions required for the existence flat bottom. The 
key his successful development was assume that the bottom amplitude 
was higher order the perturbation parameter and that the speed the 
free surface wave and the speed the free-stream flow were initially equal. 

The major limitations each these models for measuring the RST 
estuary with oscillatory flow that only steady unidirectional flows are 
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treated and that the benthic boundaries are either impermeable their motion 
must determined empirically. The solution for the motion irregular, 
erodible, and permeable bottom under oscillatory flows does not presently exist. 

Lee (13) developed stochastic theory and analyzed the particle movement 
within sand dune bed under uni-directional hydrodynamic flows. Lee (13) 
assumed that the process alternate erosion and deposition sand grains 
stochastic and that the rest periods and step lengths the sand particles 
within the profile are Gamma distributed. The Gamma distribution for the rest 
periods and the step lengths two-parameter positive semidefinite distribution. 
Once the Gamma distributions are determined for measured sand wave data 
collected estuary, the mode the distributions the rest periods, 
and the step lengths L,, could used estimate the RST the estuary. 
The celerity the sand wave would then estimated the relation 

Hino (8) developed spectral theory for sand waves also under uni-directional 
flows. derived parametric equations for the wave number and frequency 
domain spectra for sand waves based dimensional analysis. The wave 
number and the wave frequency are related the celerity for any harmonic 
the sand wave profile estuary, which given 


which harmonic frequency the wave frequency spectrum; and 
harmonic wave number the wave number spectrum the sand waves. 
The value the peak frequency the sand wave sprectrum, could provide 
average estimate the RST for the sand waves and could estimated 
Eq. from the peak wave number and the celerity. The advantage the 
spectral method Hino (8) over the stochastic method Lee (13) that 
would only necessary measure spatial record the sand wave elevations 
since the peak frequency the frequency domain spectrum may computed 
from Eq. using transformation between wave frequencies and 
wave numbers. contrast, the stochastic method Lee (13) requires both 
temporal and spatial record the sand wave profile since dispersion 
transformation presently exists between the Gamma distribution for step lengths 
and rest periods. great deal difficulty and expense involved obtaining 
time-dependent record sand wave profile under field conditions since 
the profile sensor must introduced into the flow field without significantly 
disturbing the flow. Furthermore, the oscillatory motion the wave-like 
benthic interface slowly-varying function time, then the profile sensor 
must either protected against long-term mechanical and biological degradation 
repositioned exactly the same location relatively short intervals time 
over number oscillatory cycles. This second requirement for repositioning 
may involve expensive navigational positioning systems. 

Because the potential economic savings offered the spectral method 
Hino (8) over the stochastic method Lee (13) result eliminating 
the requirement measure temporal record the sand wave profile, the 
spectral method was adopted for application the Bella-Williamson dissection 
plane approach. The essential theoretical elements from the eloquent derivation 
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the Hino spectral method (8) are succinctly summarized for convenience 
demonstrating its application. 


The random, irregular shape estuarine sand wave profiles may represented 
the harmonic analysis number sinusoids varying frequencies having 
spectral representation (compare Refs. and 22). Hino (8) derived the 
parametric equation for wave number and wave frequency spectrum for 


wave numbers; and frequencies within the equilibrium subrange for sand wave 
values are shown Fig. 


sand waves 
develope 


wave number, 


wave frequency, 


FIG. Sand Wave Spectra: (a) Wave Number Spectrum and 
(b) Wave Frequency Spectrum, 


For fully-developed (or equilibrium) sand waves, the slope the sand bed 
profile may never exceed the angle repose, the sand. This implies 
that upper limit exists for the spectral amplitude values within equilibrium 
subrange which the spectral form governed predominantly the angle 
repose and some power the wave number. For lower values wave 
numbers below this equilibrium subrange, the sand waves are still growing and 
the law equilibrium subrange fails due instability mechanism between 
the sand bed and the flowing water result the continuing transfer 
energy from the hydrodynamic shear stress the deformable sand bed. 

Wave Number Spectrum.—The wave number spectrum may estimated from 
spatial record sand wave profile. From dimensional analysis, Hino 
(see Ref. Eq. 10) has determined that the spectral density function for 


equilibrium 
subrange 
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sand wave profile within the equilibrium subrange the form 


which the wave number cycles per unit length; function 
the diameter the sand particle; the vertical elevation the sand wave 
profile above the horizontal mean any horizontal coordinate spatial sand 
wave record (see Fig. 3); and the smallest value for the wave number 
the sand wave field which the interactional instability mechanisms are 
not important. Analyzing experimental data for uni-directional flow, Hino (8) 
found that approximate value for may estimated from 


which the depth water the laboratory channel. 

Wave Frequency Spectrum.—The wave frequency spectrum may estimated 
from temporal record sand wave profile. contrast the single 
power (see Eq. throughout the wave equilibrium subrange, 
Hino (8) suggests that there are two power laws for the wave frequency equilibrium 
subrange. This dual parametric law for the wave frequency equilibrium subrange 


FIG. 3.—Definition Sketch for Sinusoidal Sand Bed 


may due, part, the form the frequency dispersion transformation 
assumed Hino (8). alternate frequency dispersion transformation which 
removes this dual parametric law for the wave frequency equilibrium subrange 

Hino (8) assumes that within the equilibrium subrange, the frequency spectrum 
follows power law’’ for higher frequencies, and power 
for frequencies slightly greater than the spectral peak frequency but less than 
the power frequencies. the higher frequency region within the 


equilibrium subrange, the frequency spectrum for the sand wave profile assumed 
from dimensional analysis given 


which the wave frequency cycles per unit time; function 
which nondimensional quantity related the critical tractive force 
and and the lower and upper bounds, respectively, the higher frequency 
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region within the equilibrium subrange for sand wave frequencies. 
The wave number, and the wave frequency, sand wave are related 
the wave celerity rate translation the waveform according 


The relationship between the frequency spectrum and the wave number spectrum 
may determined equating the areas under two differential intervals, 
and the wave number and wave frequency spectral densities, respectively. 
Equating these differential areas yields 


= ) df 


which may used estimate the wave frequency spectrum without having 
measure temporal sand wave profile. the lower wave number range 
for small values the wave numbers below both the spectral peak wave number 
and the equilibrium subrange, energy transfer between the sand waves and 
the water flowing over them become predominant since the sand waves are 
still growing. Within this lower wave number range, the celerity must considered 
function the wave number; i.e. 


Within the equilibrium subrange for wave number values Eq. 


which dimensional constant with dimensions Substituting Eq. 
into Eq. and taking derivatives both sides, obtain 


The equilibrium subrange the wave number spectrum must extend wave 
numbers which are lower than for the equilibrium subrange the frequency 
spectrum (8). The equation for the frequency spectral density function for this 
relatively lower frequency region below the wave frequency equilibrium subrange 
may now obtained substituting Eqs. and into Eq. i.e. 


which the lowest value for the wave frequency the wave spectrum 
and corresponds the lowest wave number, k,. 

summary, within the equilibrium subrange for the wave number spectrum, 
Hino (8) concludes that there must two different frequency regions for the 
corresponding wave frequency spectrum. gives the following two forms 
for the wave frequency spectrum within this corresponding wave number 
equilibrium subrange: (1) The power Eq. 10: (1/2) 
The key the successful application the spectral 
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method for estimating the RST lies determining the appropriate constant 
proportionality, and the angle repose function, 

For application the Bella-Williamson RST-OCS dissection plane method, 
was proposed use the total energy content the sand wave number spectral 
density (or, equivalently, the variance the spatial profile record) and the 
wave number equilibrium subrange defined Hino (8) measures the 
RST estuary which sand waves are present. The relative magnitude 
the total energy the sand wave profile and the equilibrium subrange amplitudes 
would provide estimates the relative energy transfer occurring between the 
hydrodynamic shear stress and the resulting deformation the migrating benthic 
interface. For comparison, data were measured both pristine estuary and 
adjacent developed (dredged) estuary. The relative ease and low cost involved 
acquiring these sand wave records are especially attractive economically and 
are fundamentally compatible with the philosophy the Bella-Williamson 
RST-OCS dissection plane approach the evaluation and preparation EIS 
for dredging operations estuaries. 


Sano Wave Data ano 


Sand wave profiles were recorded analog form submerged sonic profiler 
during two-day period. rigid steel pipe was mounted vertically the side 
boat and the sonic transducer was secured the underwater end approx 
0.8 above the sand bed. The sonic data were recorded directly magnetic 
tape analog form tape recorder aboard the boat. The sonic sampling 
locations (STA are shown Fig. 

Five buoys were implanted along the channel each the sampling stations 
shown Fig. and the distances between these buoys were measured 
means range finder (transit type) shore. The buoys were moored taut 
and anchored the bottom the estuary order minimize their drift. 
Before recording the sonic data each station, the boat was driven slowly 
along the course the buoys and the time travel between the buoys was 
recorded stop watch. The speed the boat was then computed dividing 
the known distances between the buoys the recorded time required for the 
boat cover these distances each station. The speed the boat was 
approximately equal 0.5 m/sec. The speed the boat each station was 
then maintained approximately constant value 0.5 m/sec during each 
the subsequent data recording operations. During the two-day data-recording 
period, two data recordings were taken STA four STA and one 
STA These seven data measurements are summarized Table 

order correlate the situ length the sand waves with the magnetic 
tape data, the passing each buoy the sampling boat was recorded 
voice the same channel the magnetic tape the sonic data. The pulse 
rate for the sound waves from the sonic profiler was 1/30 second. This rate 
data sampling was required greater than the speed which the boat 
was moving (approx 0.5 m/sec) order insure that sand wave elevation 
data were omitted from the recordings. order eliminate extraneous vertical 
disturbances the boat and, consequently, the data records, data were 
collected only when the water surface disturbances were minimum. order 
maintain sufficient draft below the underwater sonic transducer, data were 
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collected only during the slack periods high tide. 

These sonic sand wave data were recorded analog form magnetic 
tape and were later copied onto analog strip chart for visual observation 
and calibration. The total length each run was computed both meters and 
seconds. The elapsed time between the start and the end each magnetic 
tape record gave the total length the record seconds. The total length 
the sand wave records seconds was required order digitize the sand 
wave record analog digital converter (ADC) through PDP E10 
minicomputer. The spatial resolution, Ax, each sand wave record was computed 
multiplying the tangent the sonic profiler half-beam angle (5°) the 
vertical height the sonic transducer above the sand bed (approx 0.8 m). 
The spatial resolution, Ax, computed from the aforementioned method was 
approximately equal 0.07 This spatial resolution was then compared with 
the horizontal distance between two consecutive sonic waves insure that 
was greater than the horizontal distance between two consecutive sonic waves. 


TABLE 1.—Summary Sonic Profiler Data 


Variance, Sampling 
Length, Length, location, 


(1) (2) (3) (4) (6) (8) 


This horizontal distance was computed multiplying the average speed 
the boat the rate which the sonic waves were sent. Once the spatial 
resolution each record was determined, the total number digitizing values, 
may computed from the relation 


which the total number digitizing values; and the total length 
the sand wave profile meters. The digitizing rate required for the ADC, 
At, was then computed from the relation 


which the total time the record seconds obtained from the analog 
magnetic tape. 

The digital spatial sand wave data were first linearly detrended numerical 
algorithm and then the amplitude spectra were estimated integer finite 


91.09 236 0.0113 0.30 0.011 0.075 
127.53 471 0.1348 1.04 0.008 0.075 STA 
211.34 675 0.0615 0.70 0.005 0.038 STA 
211.34 400 0.0337 0.52 0.005 0.038 STA 
211.34 402 0.0654 0.72 0.005 0.038 STA 
216.02 540 0.0427 0.58 0.005 0.038 STA 
45.81 270 0.2382 1.38 0.022 0.017 STA 
(11) 
(12) 
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Fourier transform (FFT) algorithm base two. For real-valued sequence, 
the discrete FFT pair may given 


which the real and imaginary components the complex-valued FFT 
coefficients are identified 


order for the inverse transform exact, must also require that 


Only the real part the complex-valued sequence represents the digitized 
values the sand wave profiles. The one-sided spectral density function for 
the sand wave profile may estimated from the complex-valued FFT coefficients 
according 


which the superscript asterisk denotes complex conjugate value. The 
amplitude spectra for the sand waves recorded the three stations Coos 
Bay estuary are shown Figs. 4-9. The variance, root-mean-square wave height, 
and the equal discrete wave number interval, Ak, are tabulated Table 
The variance the zero-mean spatial sand wave profile computed from 


and the from 


N/2 1/2 


The amplitude spectra the sand wave records from the South Slough estuary, 
Coos Bay, Oregon (see Figs. 4-8) demonstrate two important points. First, 
the total energy (or, equivalently, the statistical variance the spatial sequence) 
the spectra are relatively low (e.g., the minimum 0.3 for record 
No. Fig. 4). Second, most the energy contained the fundamental 
harmonic component. This indicates that the dynamic migration these sand 
waves slowly-varying, low-energy process and that the sand bed 
slowly-undulating boundary. addition, the maximum spectral amplitudes are 
found wave numbers which lie below the smallest wave number the 
equilibrium subrange given Hino (8) for uni-directional flow the laboratory. 
This indicates that the sand waves the South Slough estuary are still growing. 
estuaries low sand wave energy, catastrophic events such storms govern 
the RST and not the hydrodynamics the average flow. the South Slough 


(14) 
(17) 
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estuary data, the rms sand wave height, varied from 0.30 

contrast, sand waves the Coos Bay main channel (see Fig. have 
higher energy content (variance the spatial sequence) dispersed over 
wider band wave numbers; however the fundamental harmonic still contains 
most the total energy. This due the stronger tidal velocities probably 
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FIG. Wave Amplitude Spectrum for Record No. 0.0113 
0.011 


2 


Wave number, k(mi') 


FIG. Wave 


Amplitude Spectrum for Record No. 0.0615 
0.005 


augmented the ship propeller disturbances the main channel. Since this 
spectral comparison only visual comparison, the spectral band width parameter 
which provides estimate the width the wave frequency spectrum (compare 
Ref. 3), should computed from the wave frequency spectrum order 
obtain the dispersion the spectrum about the spectral peak. The sand 
wave frequency spectrum may obtained from the wave number spectrum 
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using theoretical empirical transformation between the wave numbers 
and the wave frequencies. flow chart for obtaining the RST from wave 
number spectrum using dispersion transformation shown Fig. 10. 

For the sand wave record from the Coos Bay main channel (see record No. 
Fig. 9), the wave number the maximum spectral amplitude falls very 
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> 


Wave number, k (mi!) 


FIG. Wave Amplitude Spectrum for Record No. 0.0337 
0.005 
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-050 
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FIG. Wave Amplitude Spectrum for Record No. 0.0654 
0.005 


close the minimum wave number value for the equilibrium subrange given 
Hino (8). note that the minimum wave number value for the equilibrium 
subrange, given Hino (8) slightly less than the first harmonic component. 
This demonstrates that some additional care must taken when digitizing analog 
data order insure that the value for the minimum wave number for the 


equilibrium subrange, will much greater than the fundamental harmonic 
component. 
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For example, assume that will appear.in the mth harmonic (e.g., 


Ak; 1), then the equal discrete wave number interval, Ak, may 
estimated from Eq. according 


From this value and the spatial resolution, Ax, computed from analysis 


Spatial Spectrum, lal? (k) (m? x 10) 


k, = 0.0375 mi! 


Wave number, k (mi!) 


FIG. 8.—Sand Wave Amplitude Spectrum for Record No. 0.0427 m?; 
0.005 


Swi) k 


Wave number, k (mi!) 


FIG. Wave Amplitude Spectrum for Record No. 0.2382 
0.022 


the data acquisition system, the total number digitized values, may 
computed from: 


For the Coos Bay main channel, assume that will appear the fifth 
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harmonic, then the equal discrete wave number interval, Ak, may estimated 
from Eq. 0.0033 Substituting this equal discrete wave 
number interval into Eq. yields approx 4,096. Therefore, the spatial 
sand wave data should have been digitized with least 4,096 values instead 
the 1,024 values used the present analysis, order insure that 
would much greater than the fundamental harmonic component and better 
resolve the equilibrium subrange. 


FIG. 10.—Flow Chart for Estimating RST from Wave Number Spectrum Dispersion 
Transformation 


These amplitude spectra demonstrate the significant differences between the 
hydrodynamic processes pristine estuary (South Slough) and developed 
estuary (Coos Bay Main Channel). The RST the pristine estuary low-energy 
process, which storm driven, while the RST the developed estuary 
high-energy process driven the hydrodynamics the average flows. 
Additional laboratory and field measurements are required, however, verify 
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the constants proportionality for the equilibrium subrange and the dispersion 
transformation. 


The equilibrium subrange and the total energy (variance the spatial sequence) 
the wave number spectrum for sand waves estuary may used 
quickly and economically locate the RST the Bella-Williamson dissection 
plane method for diagnosing the chronic environmental impacts dredging 
estuaries. spatial record sand wave profile estuary may 
easily measured underwater sonic profilers and the wave number spectral 
amplitude may quickly estimated FFT algorithm. comparison between 
the rms wave height (which may determined from the total energy content 
computed the FFT coefficients) with the equilibrium subrange the theoretical 
wave number spectrum developed Hino (8) provides reasonable estimate 
for measure the RST. Some care must taken both recording and 
digitizing spatial sand wave data and the important analytical relationships between 
the spatial resolution the data acquisition equipment, the speed the sampling 
boat, and the minimum value the equilibrium subrange wave number are 
presented. 

The method described the present analysis may extended compute 
the wave frequency spectrum without requiring the measurement temporal 
sand wave profile provided that dispersion transformation between wave 
frequencies and wave numbers known. uniform power law’’ wave 
frequency spectrum which follows the assumption that the wave 
celerity proportional the bed shear velocity. Additional analytical, 
experimental verifications, both, are required determine the dimensionless 
constant oscillatory flows. 

The potential economic savings which may realized the preparation 
and evaluation EIS the judicious application the Bella-Williamson 
dissection approach are tremendously encouraging. key ingredient their 
approach the successful integration disciplinary data and analysis which 
may rapidly and economically collected and analyzed. The spectral analysis 
sand wave profiles FFT estuary presented within the context 
disciplinary input the Bella-Williamson dissection plane method. Details 
regarding the integration disciplinary data are given Bella and Williamson 
(2). addition, recommendations are presented for extending the present 
knowledge the wave number equilibrium subrange and the wave frequency 
dispersion transformation for sand waves under oscillatory flows. 


ACKNOWLEDGMENTS 


are especially grateful for the unique opportunity provide disciplinary 
data Bella and Williamson for integration the development 
their pioneering work the RST-OCS dissection plane approach for diagnosing 
the chronic impacts dredging estuaries which was sponsored the National 
Science Foundation-Research Applied National Needs (NSF-RANN) under 
Grant No. ENV71-01908. Equally rewarding and professionally stimulating have 
been the collaborative opportunities with the following NSF-RANN team 


FEBRUARY 1981 HY2 


Stander. The sonic profiler used the data acquisition was made available 
Nordin the United States Geological Survey, Department the Interior. 
Additional financial support (for RTH) was provided the Oregon State 
University Sea Grant College Program, National Oceanic and Atmospheric 
Administration, Office Sea Grant, Department Commerce, under Grant 
No. 04-6-158-44004. dedicate this effort memory our late NSF-RANN 
team member McCauley, benthic biologis:, educator, and good man. 
also pleasure extend our appreciation the reviewers for their valuable 
suggestions how improve both the clarity and accuracy our efforts. 


The equilibrium wave frequency spectra derived Hino (8) follows 
power law’’ for the higher frequencies within the equilibrium subrange and 
power for frequencies slightly greater than the spectral peak frequency 
but less than the power frequencies. This indicates that two frequency 
dispersion relationships are required for the sand waves within the wave number 
equilibrium subrange. Since the wave frequency spectrum prime importance 
the estimation the RST, uniform ‘‘—3 power developed for 
the entire wave frequency equilibrium subrange. This results one-to-one 
correspondence with the wave number equilibrium subrange. 

Following the derivation for the equilibrium subrange for surface gravity waves 
Phillips (16), note two important points. First, the constants proportiona- 
lity used his dimensional analysis are dimensionless; and, second, the same 
dimensional constant used the parametric equation for the wave frequency 
spectrum also appears the frequency dispersion equation (i.e., the gravitational 
constant). shall incorporate these two points our following formulation. 

Gradowczyk (6) and Hino (see Ref. Eq. 16) have shown that the celerity 
the sand wave profile proportional the bed shear velocity. We, therefore, 
assume the following for the frequency dispersion transformation given Eq. 


which dimensionless constant proportionality. Eq. now becomes 


Substituting this differential relationship into Eq. yields 


The lowest frequency, now corresponds the lowest wave number, k,, 
the equilibrium subrange given Hino (8). 

The upper wave number range the equilibrium subrange assumed 
Hino (8) given 


which the diameter the sand particle. Accordingly, the upper frequency 
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and uniform for the entire wave frequency equilibrium subrange 
now given 
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The following symbols are used this paper: 


complex finite Fourier transform (FFT) coefficient; 

celerity sand wave; 

diameter sand particle; 

sand wave frequency; 

lower and upper frequency bounds for the power 
equilibrium subrange (8); 

lowest frequency the equilibrium subrange; 

peak frequency the wave frequency spectrum; 

function 

root mean square wave height; 

depth water channel; 

imaginary part complex FFT coefficient; 

imaginary unit number; 

wave number 

smallest wave number the equilibrium subrange; 

peak wave number the wave number spectrum; 

length sand wave record; 

wave length sand wave; 

total number digitizing values; 

real part complex FFT coefficient; 

spectral density function for the sand wave profile; 

total time the sand wave record; 

wave period sand wave; 

variance the sand wave spatial profile; 

shear velocity; 

horizontal coordinate axis; 

instantaneous vertical elevation sand wave profile measured 
positive upwards from horizontal bottom; 

function 

dimensionless constant proportionality for frequency disper- 
sion; 

dimensional constant proportionality for frequency dispersion 
(8); 

equal discrete wave number interval; 

equal discrete time interval; 

equal discrete spatial interval; 

fluid density; 

density sand particles; 

angle repose sand particle; and 

nondimensional quantity related the critical tractive force 
the sand bed. 
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and Ramesh 


Many practical problems connected with alluvial rivers and canals require 
knowledge the sediment load carried them the load they can safely 
carry without danger aggradation degradation. For example, the loss 
storage capacity reservoir built alluvial river primarily dependent 
the sediment load carried the river and transport capacity flow 
various sections the back water reach. Also, the design nonscouring, 
nonsilting alluvial canal requires that only much sediment allowed into 
the canal can carry under the given hydraulic conditions. Essentially, 
therefore, the hydraulic engineer required predict the sediment load carried 
channel under equilibrium conditions for given sediment, fluid, and flow 
characteristics. 

Wash load and bed-material load are the two components the sediment 
load transported alluvial channel (5). known that the wash load 
poorly related the hydraulic conditions the stream, and few attempts 
have been made the past predicting this part the load, despite the 
fact that may form major part the load many channels. The bed-material 
load function the sediment, fluid, and flow characteristics and 
component amenable rational analysis. The bed-material load may move entirely 
bed load close the bed, part may move suspension suspended 
load, the bed material fine and the shear stress large enough. Hence the 
total load implied this paper the sum the bed load and suspended 
load those sediment fractions originating the bed; other words, does 
not include the wash load. 

Total load computations may made (6) either computing the bed load 
and suspended load separately and adding them up, using approach 
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TABLE 1.—Data Used Analysis 


Median | 


Range 
size of Numbe 


chosen per secon percentage 


Source and Maxi- Maxi- 
reference random mum mum 


(2) (5) (11) 


Kalinske and Hsia (9) 0.100 
Laursen (9) 0.117 
Vanoni and Brooks (9) 0.330 
Nomicos (9) 0.250 
Laursen (9) 0.210 
Sato and Kikkawa (9) 0.0121 
(11) 0.100 
Znamenlkay (9) 0.163 
Nomicos (9) 0.390 
Vanoni and Brooks (9) 0.278 
Brooks (9) 0.310 
Nomicos (9) 0.340 
Southampton, A (9) 0.251 
Barton-Lin (9) 0.170 
Straub (9) 0.235 
USGS (9) 0.092 
Vanoni-Hwang (9) 1.200 
Znamenskay (9) 0.205 
Vanoni-Hwang (9) 0.248 
Franco (9) 0.169 
Mutter (9) 0.360 
USGS (9) 

USGS (9) 

Gilbert (9) 

USGS (9) 

O’Brien (9) 

UBC (9) 

Stein (9) 

USWES (9) 

USGS (9) 

USGS (9) 
Southampton, B (9) 
Gilbert (9) 

USGS (9) 

Bharat Singh (9) 
Gilbert (9) 

USGS (9) 

Williams (9) 

Liu (9) 

Liu (9) 

Casey (9) 

Liu (9) 

(9) 
Gilbert (9) 

Liu (9) 

Liu (9) 

Corps Engineers (9) 


Colorado River 
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TABLE 1.—Continued 


Niobrara River (3) 
Leopold (9) 

Mountain Creek (9) 
Portugal River (9) 
Middle Loup (9) 


and canal data. 


which the total load directly related the hydraulic and sediment 
characteristics. review the existing methods both these categories 
available (6) and suggests the need for additional work this field view 
the unsatisfactory accuracy these computational methods. Hence the writers 
undertook investigation with the dual objectives testing the accuracy 
some the existing computational methods and proposing alternative 
and, possible, improved computational technique. The results that 
investigation are presented herein. 


Data 


Since large amount data collected flumes and natural streams under 
equilibrium conditions was available, was felt unnecessary carry out any 


experiments. These data from alluvial rivers, canals, and flume studies were 
carefully selected include only those which information concerning the 
wash load was available, thereby assuring that the total load used the analysis 
excluded the wash load. The data cover wide range sediment size and 
other hydraulic conditions; only data from ripple-dune and transition regimes 
were chosen for analysis. Antidune data were not included because this regime 
occurs seldom the field. The sources these data and some other relevant 
details are given Table Natural sediment with relative density 2.65 
was used all these studies. 

correction for bank friction has been applied the river data, and the 
depth flow has been taken the hydraulic radius with respect 
the bed R,. The bed and sides canals have been assumed the 
same nature, and the hydraulic radius was used for R,. The value 
was computed for the flume data adopting the procedure side wall correction 
given Einstein (4). 


Concert Errective Stress 


Vittal, al. (10) found that the dimensionless sediment transport parameter 
was closely related the dimensionless shear stress for flow with plane 
bed, i.e., for flow which the bed free from undulations. The relationship 
obtained them shown Fig. and 
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which rate total load transport (by weight) per unit width; 
gravitational acceleration; median size the sediment; mass density 
the sediment and the fluid, respectively; and average shear stress 
the bed given y,R,S which hydraulic radius with respect the 
bed, and slope. worth mentioning that the data used developing 
Fig. showed very little scatter this plot (see Ref. 10). Hence the writers 
felt that Fig. provides good basis for further analysis and adopted this 
further analysis. The rather high dependence the transport rate the 
shear stress low values the shear stress (see Ref. for instance) typical 
most transport relations based the shear stress. Fortunately, however, 
not range great practical interest. 

Vittal, al. defined the shear stress for ripple and dune beds the 
shear stress required give the same total load transport the same sized 
material plane bed. Inasmuch some the shear stress undulated 


10 


0-01 
0-00002 0.0001 0.00! 0.01 


FIG. 1.—Vittal, al. Sediment Transport Relation for Plane Bed (10) 


bed ineffective causing sediment motion, the effective shear stress always 
less than the average shear stress the bed. accordance with the foregoing 
definition, the effective shear stress can obtained from Fig. for the known 
value sediment transport rate (or undulated bed. Analysis 
data from ripple and dune bed regime indicated (10) that uniquely 
The hydraulic radius corresponding grain roughness given (in units) 


which the Stricker’s roughness coefficient, and mean velocity flow. 
The relation between and was proposed used for determination 
and subsequent determination from Fig. Vittal, al. (10). 
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During the course this investigation check was made the adequacy 
the and relation. The data listed Table indicated that 
there unique relation between the foregoing parameters, and accuracy 
predicting Fig. significantly affected the errors the prediction 
Therefore the writers attempted provide better predictor for 
the following lines. 

Functional Relationship for al. (10) have shown that greater 
than case ripple and dune beds, suspended load present. This 
was ascribed the fact that part the form drag would effective throwing 
material into suspension, although would not help bed-load movement. 
and fall velocity the sediment still fluid), that is, when there 
suspended load movement and yet the bed shear stress greater than the 
critical shear stress leading bed-load movement, would equal 
when the bed plane, tends and tends but when undulations 
are present may expected depend also. Thus the functional 
relationship based this physical reasoning may written 


o o w 


Eq. may also obtained from principles dimensional analysis follows: 
ing that depends p,, and and also that related through 
Chezy-type resistance equation, one may write which 
obviously reduces Eq. 

Variation with and fall velocity the sediment 
has been calculated using Rubey’s equation, that 


The data listed Table were then grouped into different ranges u,/w 
and Fig. was prepared accordance with Eq. can seen clearly that 
increases with increase u,/w for any given This rather 
significant because generally increases with increase and the increase 
towards suspension increases. Further, decreases with increase 
for any value Detailed examination Fig. indicates that the relationship 
for can expressed 


equal unity; other words known (6) that suspended load 
negligible when 0.5. Thus 
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The values determined from Fig. for other u,/w values were plotted 
against (not shown here) and the relationship obtained 


for 0.5 
WwW 

The values were calculated for the data using Eqs. and 
computed. The agreement the data with Eq. shown Fig. The over-all 
error using Eqs. and well Fig. for the prediction the transport 
rate apparent from Fig. (to considered subsequently herein). 

Eq. can also written 


Since function one can write 


Eq. valid for ripple-dune well plane bed regimes; the latter case 
and are equal and the term within the large bracket Eq. merely 
reduces t,. Note that the same term reduces when suspended load 
absent, since such case. The data listed Table along with 
the plane bed data used Vittal, al. (10) are plotted Fig. accordance 
with Eq. 11. The values were calculated from Eqs. and The good 
correlation between the two parameters Fig. may expressed 


the range 


The range given Eq. should cover most practical situations. One can 
use the mean line drawn Fig. outside the foregoing range 
along with Eqs. and affords simple means calculating the total 
load. The error computation this procedure was found within 
+40% for over 80% the data. 


The various sets data listed Table were used check the accuracy 
few the many relations for transport proposed recent years (1,2,3,5,7,8, 
12,13). not the intention herein check the accuracy all the existing 
transport relationships. fact, hoped that such checks some relationships 
using the data listed this paper well check the writers’ relation 
using additional data will carried out the discussers this paper. Some 
conclusions could also drawn comparison one relation with another. 
Hence only the following relations were chosen for such check: (1) Ackers 
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and White (1); (2) Engelund and Hansen (5); and (3) Graf and Acaroglu (7). 
Ackers-White Relation.—The relationship for load transport proposed 
Ackers and White can written 


Ys d C; 
inwhich 
Ay, 10D 


and total load concentration weight. 
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FIG. 5.—Observed Transport Rate Plotted Against Transport Rate Calculated 
Ackers-White Method (1) 
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0.23 


3 


For 60.0 
C,=0, C,=0.025, C,=0.17 and C,=1.50 


which being the kinematic viscosity. The values 
were calculated for all the data using Eqs. through 17, and hence 
was computed. The computed values are plotted against the observed trans- 
port rates Fig. The error the predicted transport rate within +50% 
for 80% the data. 


Engelund-Hansen and Hansen gave the total load equation 

Fig. shows plot against for all dune-bed data. may seen 


that, general, the sediment transport rate much greater than that given 
Eq. 18, and that 


fo,=0.25 


Graf-Acaroglu Relation.—The sediment transport law Graf and Acaroglu 


a 


This can expressed 


Garde (6) found that large amount flume and field data followed the relation 

Comparison Eqs. and shows that possibly different sets data have 

been used the development these equations, and neither them satisfies 


all the data. One might also argue that the sediment transport function more 
complex than indicated either Eq. Eq. 22. 


The total load transported alluvial channel can expressed functional 
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The following symbols are used this paper: 


constants; 


total load concentration weight; 
depth flow; 

median diameter sediment; 
Darcy-Weisbach friction factor; 

gravitation acceleration; 

constant; 

Manning’s roughness coefficient; 
total sediment rate, weight per unit width channel; 


3 


Available data from canals, rivers, and flume studies, for regimes other than 
the antidune regime, indicate that Eq. 12, which based the foregoing 
functional relationship, satisfactory sediment transport relation for the 
commonly encountered regimes flow. This relationship appears more 
accurate than some the existing total load equations. 


ALLUVIAL CHANNELS 


hydraulic radius with respect grain roughness; 
energy slope; 

average velocity over the cross section; 

shear velocity; 

unit weight fluid; 

unit weight sediment; 

dynamic viscosity; 

kinematic viscosity; 

mass density fluid; 

average bed shear stress; 

average bed shear stress corresponding grain roughness; 
effective shear stress; 

dimensionless grain shear stress; 

dimensionless total load transport; and 

fall velocity sediment still fluid. 
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INTRODUCTION 


The serious consequences failure dam structure make imperative 
that these structures designed appropriately for all possible hazards. The 
occurrence strong earthquakes and the resulting damage that had occurred 
some structures has drawn attention hazards posed earthquakes. 

Construction the concrete-arch Green Lake Dam one North America’s 
highly seismic areas required assessment the hazards and risks associated 
with this level seismicity. This was evaluated determination anticipated 
bedrock acceleration the dam site, determination the response behavior 
the dam, and detailed stress analysis ensure that the stress levels 
the dam would within generally acceptable limits. 

The Green Lake Hydroelectric Project located Southeast Alaska about 
air miles (19 km) southeast Sitka the west side Baranof Island 
(see Fig. 1). The project being constructed the city and borough Sitka, 
Alaska. The arch dam located across the Vodopad River the mouth 
Green Lake, existing glacial lake approximate El. 240 (73 m), which 
drains tidewater Silver Bay some 2,000 (610 downstream. The reservoir 
will have total storage capacity 90,000 acre-ft (111 10° the normal 
maximum reservoir, El. 390 (119 m). Water for power generation will 
conveyed the powerhouse tidewater 8-ft (2.4-m) diam concrete-lined 
tunnel, 1,900 (580 long. The powerhouse will contain two vertical axis 
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Francis turbines connected two generators with total rated capacity 
16,500 kW. 

The dam, currently under construction, will concrete double curvature 
arch, shown Fig. The crest length along the axis will 450 (137 
m), subtending angle 110°. The dam will 210 (64 high, 
(1.8 thick the crest, and (4.9 thick the base the crown 
cantilever. The total volume concrete the dam will 26,000 (19,900 

Seismic records indicate that approx earthquakes have been recorded within 
100 miles (160 km) the site between 1899 and 1976. The records suggest 
that least one those earthquakes may have caused significant ground motion 
the dam site. Recognizing that the record seismicity Southeast Alaska 
very short, and likely deficient with some inaccuracies, was concluded 
that study was necessary realistically assess the seismic hazards and risk 
the site and that the design the arch dam should take into account the 
anticipated seismic ground motions identified the study. 


Green Lake dam site located the west side Baranof Island. Foundation 
rock for the dam very competent, slightly metamorphosed greywacke. The 
topography precipitous, rising approx 4,000 (1,200 within 1-1/2 miles 
(2.4 km) the shoreline. The core the island consists largely Jurassic 
and Cretaceous-age, before present (ybp), granitic 
rocks which have intruded and metamorphosed older stratified sedimentary rocks 
the east and west sides the island. Subsequent deformation the rocks 
occurred, including folding and movement strike slip and vertical motion 
along faults. The period major deformations probably terminated more than 
10,000,000 ybp. 

Baranof Island was intensely glaciated during Pleistocene time (between 10,000 
ybp). The continental ice sheet covered the island about 
3,000 (915 m). The last ice probably disappeared about 10,000 ybp leaving elevated 
terraces, rounded rock lips (Green Lake dam site), U-shaped valleys (Green 
Lake drainage basin), and deep V-shaped fiords (Silver Bay). 


AND 


Southeast Alaska presently the site major plate convergence and its 
associated tectonics. The plate boundary between the North American plate 
and the Pacific plate along the Fairweather-Queen Charlotte Islands Fault 
which lies about miles (56 km) southwest the dam site (see Fig. 3). That 
fault has been associated with numerous earthquakes, several which fall 
into the (8+ magnitude) category. 

Baranof Island bounded three major faults, the Fairweather-Queen 
Charlotte Islands the southwest, the Peril Strait the northeast, and the 
Chatham Strait-Lynn Canal the east and southeast (see Fig. 3). The Chichagof- 
Sitka Fault traverses part the island near Sitka and its trace passes within 
1/2 mile (0.8 km) the dam site. These faults were considered the 
major tectonic along which, earthquakes occurred, potentially 
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damaging ground motion could occur the dam site. Therefore, they were 
the subject evaluation during the study. For purposes this study fault 
was assumed capable there was evidence movement within the last 
35,000 yr. 

Clearly, the Fairweather-Queen Charlotte Islands Fault the most important 
capable fault Alaska since acts the boundary between the North American 


and Pacific plates. The fault assumed long 470 miles (752 km) 
(6,20). 


EXPLANATION 


SYMBOL MAGNITUDE 
7.48.5 


FIG. 3.—Earthquake Epicenters and Faults 


The Peril Strait Fault about 110 miles (176 km) length (6). Major movement 
the fault appears have occurred prior the Miocene age (25,000,000 
ybp) (13); however some low level seismicity recorded the north could 
associated with the fault and therefore was assumed capable. 

The Chatham Strait-Lynn Canal Fault was conservatively classified capable 
circumstantial evidence consisting small amount microseismicity which 
could associated with the fault. Also, toward the south end, the literature 
reports interpretations marine geophysical surveys which indicate deformations 
and faulted sediments the Holocene age (less than 10,000 ybp) along its 


138° 136° 134° 132° 
| 
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trace (26). The fault believed some the southeast continuation 
the major active Denali fault system the north, giving combined length 
greater than 265 miles (424 km) (13). 

The Chichagof-Sitka Fault has three major segments with combined length 
between 100 miles (160 km) and 125 miles (200 km) (13). The major movement 
occurred over 25,000,000 ybp. There known geologic evidence historic 
seismicity suggest activity during the past 35,000 yr. Therefore, this fault 
was considered noncapable. short fault which splays from the Chichagof-Sitka 
Fault and passes through Green Lake about 1/2 mile from the dam site was 
carefully inspected from the air and the ground. Evidence indicate movement 
within the past 35,000 along that fault other possible capable faults was 
not detected; thus they were considered noncapable. 

Evaluation fault capability included aerial flights different sun angles 
for evaluation lineaments and detect there were other significant faults 
which had geomorphic evidence young movement. Geologic mapping and 
foundation exploration the dam site, well the aerial reconnaissance, 
revealed faults passing through the foundation; therefore permanent dam 


foundation ground displacement was not seismic consideration for the design 
the arch dam. 


Maximum 


The maximum credible earthquake (MCE) for given fault defined 
the largest magnitude event, independent time, which can reasonably 


postulated occur based upon existing geologic and seismologic evidence. 
Lacking other parameters, maximum credible earthquakes were estimated for 
the fault rupture length assuming that one half the surface mapped length 
could rupture during single event. Based the maximum credible earthquake 
determined for each capable fault, the maximum horizontal ground acceleration 
the dam site was established using the attenuation curves Schnabel and 
Seed (22); duration shaking was determined from Bolt al. (2) and Housner 
(10). summary the maximum credible earthquakes and estimated rock 
motion parameters for the dam site given Table 


TABLE 1.—Green Lake Dam Site Maximum Credible Earthquakes and Horizontal 
Rock Motion Parameters 


Rock Surface Motion Parameters 


Minimum site rock Duration 
distance, earthquake acceleration strong shaking 
Fault name miles magnitude (22), seconds (2,10) 
(1) (4) (5) 
Chatham Strait- 
Lynn Canal 
Peril Strait 


Fairweather-Queen 
Charlotte 
Islands 


q 
(29) 
0.25 
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Although the frequency earthquakes from the Fairweather-Queen Charlotte 
Fault greater, MCE from the Chatham Strait-Lynn Canal Fault would cause 
the maximum peak bedrock acceleration the dam site. 


Risk 


Seismic risk evaluations were carried out after the methods outlined Marachi 
and Dixon (14). Earthquakes occurring within 100-mile (160-km) radius circle 
about the site (see Fig. 3), formed the data bank for the historic analysis. 
Two methods analysis proved valuable comparing the seismic risk terms 
return intervals; those methods were the site region model and fault model. 
The site region model assumes the random occurrence earthquakes space 
and the independence events; results from this method should regarded 
average for particular area. Attenuation areas were developed based 
the curves Schnabel and Seed (22). 


+ 


0.0 0.2 0.3 0.6 
MAXIMUM ROCK ACCELERATION, 


AVERAGE RETURN PERIOD FOR EXCEEDING 


FIG. 4.—Average Recurrence Horizontal Rock Surface Acceleration 


The fault model only included data from the Fairweather-Queen Charlotte 
Islands Fault Zone due the lack sufficient geologic/seismic information 
for other regional faults. The magnitude recurrence relationship was based 
earthquakes instrumentally recorded along the 470-mile (752-km) length the 
fault. 

Recurrence intervals for the two methods analyses were interpreted (see 
Fig. 4). The site region model assumes the random distribution earthquakes 
both time and space and best suited analyze seismic environment 
which earthquakes are uniformly distributed over relatively large areas. 
Conversely, the fault model assumes that the earthquakes originate only 
known fault traces. Fig. indicates that recorded earthquakes are not uniformly 
distributed throughout Southeast Alaska and there preference for most 
earthquakes associated with known faults. The fault model appeared best 
suited for evaluation the seismic risks for the Green Lake Dam, even though 
recognized that epicenter locations could mislocated along single fault 
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trace and the Fairweather-Queen Charlotte Islands Fault Zone could miles 
(32 km) more wide adjacent the site. 


Two earthquakes were selected for consideration design. They are termed 
the Class and Class earthquakes. The Class earthquake the MCE. 
The dam was designed ensure that the structure could withstand the loads 


TABLE 2.—Green Lake Dam Site Horizontal Bedrock Ground Motion Parameters 


Earthquake parameter 
CLASS 
Fault Name 


Fairweather-Queen 
Variable Chatham Strait Charlotte Islands 
(1) (2) (3) 
Richter magnitude 
Site distance fault, miles 
(kilometers) 
(53) 
Maximum peak rock acceleration, 0.23-0.30 
Probability exceedance 100 yr, 
percentage 
Duration strong ground shaking, 
seconds 


GROUND ACCELERATION 
(FT PER 


TIME (SECONDS) 


FIG. 5.—Class Earthquake Time History 


from the Class earthquake without uncontrolled release water from 
the reservoir. The Class earthquake can reasonably expected occur 
once the average every 100 during the life the project. The dam was 
designed withstand the loads from the Class earthquake and remain 
functional. 

The important characteristics ground motion are the amplitude accelera- 
tion, duration strong shaking, and the frequency with which the acceleration 
traces cross the time axis which may cause stress reversals. Table summarizes 
the important characteristics the Class and Class earthquakes selected. 
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The frequency with which the direction acceleration changes may seen 
Fig. From these rock motion parameters, accelerograms were developed 
based the simulated M8+ ground motion accelerogram reported 
Jennings (11). The Class and Class ground accelerations were scaled 
0.4 and 0.23 respectively. Both upstream-downstream and cross-canyon 
horizontal components were considered. The time history for one component 


horizontal rock surface motion the Class earthquake presented 
Fig. 


Hydrodynamic loads were found significant portion the total dynamic 
loads. The method used account for hydrodynamic interaction the widely 
accepted approach developed Westergaard. Recent studies 
model arch dam subjected low-level excitations have shown reasonable 
agreement between measured pressures and those obtained using Westergaard’s 
solution (17). The volume water considered interact with the dam had 
parabolic cross section with base width equal 7/8 the height the 


FIG. 6.—Hydrodynamic Forces: (a) Upstream-Downstream Motion; and (b) Cross- 
Canyon Motion 


dam. Hydrodynamic forces were applied the upstream-downstream and 
cross-canyon directions (see Fig. 6). 


The damping structure vibratory motion difficult evaluate and 
inclusion damping effects dynamic analyses must based experience 
and judgement. Shaking tests evaluate damping have been performed for 
concrete dams Japan, Italy, and the United States, using both model and 
prototype structures. Low level excitations produced small displacements. 
Damping ratios from 6%-10% were obtained from tests arch dams using 
large excitations (4,7). These tests indicate that energy losses increase 
displacements become large. Because the anticipated opening contraction 
joints the Green Lake Dam during vibration, damping ratio 10% was 


judged appropriate for the magnitude the ground motions for which 
the dam was designed. 


AND Properties 


Modeling foundation-arch dam interaction very complex because 


» 


h! 
(a) 
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arch dam’s three-dimensional geometry. Studies made for Auburn Dam showed 
that rock modulus had little effect the dynamic characteristics the dam 
(25). Green Lake foundation hard and not complex and infinite 
instantaneous rock modulus was considered appropriate for the dynamic analysis. 

The increase concrete modulus during dynamic loading well documented 
laboratory tests. instantaneous concrete modulus elasticity 4,000,000 
psi (27,560 MPa) was used, compared the static modulus 3,000,000 
psi (20,670 MPa). Dynamic strength values are dependent the rate loading. 
Under dynamic testing, compressive strengths 20%-30% above static values 
and tensile strengths 30%-66% over static values have been cited (25). 


comparison between static and dynamic values used the design shown 
Table 


TABLE 3.—Green Lake Dam Static and Dynamic Properties 


Variable, pounds per square inch Dynamic 
(1) (3) 
Rock deformation modulus Infinite 
Concrete modulus elasticity 4,000,000 
Concrete ultimate compressive strength 4,800 
Concrete ultimate flexural strength 540 


Typical cantilever element 
Crown cantilever 


arch element 
CANTILEVER NO. 


405 
360 
320 
280 


240 
215 


FIG. 7.—Model Dam 


Horizontal construction joints were considered remain uncracked provided 
that the ultimate flexural strength was not reached, but vertical contraction 
joints were assumed have negligible tensile capacity. 


The trial load method was utilized the design analysis. The dam was divided 
into system cantilever and arch elements (see Fig. 7), each set elements 
occupying the entire volume the dam. During vibration, the dam was considered 
deflect dynamically group fixed end arches and vertical cantilevers. 
The loads the dam are divided between the arch and cantilever elements 
such manner that the strains the common nodes are identical. 
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The dam was designed for the following conditions: (1) normal full reservoir 
with low concrete temperatures; (2) minimum design reservoir elevation with 
high concrete temperatures; (3) absolute minimum reservoir with 
average concrete temperatures; (4) probable maximum flood with average 


concrete temperatures; (5) design earthquake with Case and (6) maximum 
credible earthquake with Case 


Anatysis 


general, for structures with uncomplicated geometry, the higher vibration 
modes contribute less the total response than the lower modes. modal 
sensitivity analysis was performed determine the lowest mode vibration 
above which significant change stress encountered. Table the 
relation stress value the number modes considered shown for 
node common Arch and Cantilever The stress based 12-mode 
analysis was normalized seen that significant stress change occurred 


TABLE 4.—Green Lake Dam Effect Number Modes Stress Values 


Number modes 


(1) 


Stress value 


(2) 


when modes order higher than were included. Consequently, the design 
analysis was limited modes. 


time history analysis arch dam results enormous amount 
data which have evaluated systematically. particular interest are: (1) 
The maximum compressive and tensile stresses; (2) whether high stresses occur 
over large regions isolated areas; (3) the duration stress above certain 
level; (4) whether the tensile stresses occur the upstream downstream 
face; and (5) the magnitude and nature the other stresses occurring the 
same region the same instant. 

The first step the evaluation was list the compressive and tensile stresses 
descending order, and also their locations and times occurrence. The 
next step was plot stress-time histories points where stresses were significant. 
From these plots the frequency and duration above any given stress level were 
determined. critical instants, stress were taken. These 
showed the stresses over the entire dam given instant and areas 
high compressive and tensile stresses were mapped. Also, taking stress 
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successive instants, the pattern stress change over the entire 
dam was obtained. Adjustments the dam configuration were made until the 
design analysis results indicated acceptable performance. 


Stresses Due Crass 


Since ground motion varies function time, stresses various points 
the dam also vary with time. Table shows the arch and cantilever stresses 


TABLE 5.—Class Earthquake Stresses Time When Maximum Compression Occurs, 
+Compression, —Tension, pounds per square inch 


Arch Elevation, Cantilever Number 
number feet 


(1) (2) 


Upstream 788 
Downstream 
Upstream 923 (Maximum) 
Downstream 
Upstream 866 
Downstream 
Upstream 769 
Downstream 
Upstream 550 
Downstream 


Upstream 

Downstream 

Upstream 

Downstream 

Upstream 

Downstream 

Upstream 

Downstream 

Upstream 106 
Downstream 
Upstream 
Downstream 354 (Maximum) 


the upstream and downstream face the dam the instant when maximum 
compression occurred. Similarly, Table shows the stresses when maximum 
tension occurred. The maximum compressive stress 923 psi (6.36 MPa) (arch 
stress), shown Table and the cantilever tensile stress 132 psi (0.91 
MPa) Table indicate that adequate factor safety was attained. The 
computed arch tensile stresses shown the tables need further explanation 
for proper evaluation. 

Since vertical contraction joints have negligible tensile capacity, the computed 


(a) Arch Stresses 
405 254 173 
166 187 
360 463 207 
328 241 
320 461 283 
335 102 
280 264 
350 
240 
269 
Cantilever Stresses 
405 
—39 
320 117 237 
280 
240 
181 
195 
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negative arch stresses merely indicate that the particular vertical contraction 
joints were partially open the instant reflected the analysis. Physically, 
when the joints partially open for very short periods (fractions seconds) 
the stiffness the arch somewhat reduced and the arch load redistributed 
cantilever load for fraction second. Fortunately, this condition only 
occurs the upper arches where the cantilever stresses are small and therefore 
the redistributed arch load does not result excessive cantilever stresses. Based 
the results shown Tables and was concluded that the dam would 


TABLE Earthquake Stresses Time When Maximum Tension Occurs, 
+Compression, —Tension, pounds per square inch 


Elevation, 
in feet 


Upstream 
Downstream 


—358 (Maximum) 


Upstream 678 
Downstream —204 
Upstream 635 
Downstream 
Upstream 582 
Downstream 
Upstream 465 
Downstream 


Upstream 


Downstream 0 
Upstream 210 
Downstream —132 (Maximum) 
Upstream 186 
Downstream -71 


Upstream 
Downstream 
Upstream 
Downstream 
Upstream 
Downstream 


not sustain damage under the Class earthquake which would disrupt normal 
operation the project. 


When considering the effects the loading 0.40 Class earthquake, 
should borne mind that its occurrence would extremely rare 
event and that minor damage the structure can tolerated. However, the 
damage would not result uncontrolled release water from the reservoir. 
The maximum compressive stress for the Class earthquake about 40% higher 
than for the Class II, and the stress pattern similar. Fig. shows the arch 


Cantilever Number 
number 
(1) (7) (8) 
(a) Arch Stresses 
360 245 171 
223 108 
320 345 179 
222 
280 249 
241 
240 
221 
; (b) Cantilever Stresses 
405 
360 
320 
233 160 
240 175 165 
129 134 
195 
303 


ARCH DAM 


TABLE Earthquake Stresses Time When Maximum Compression Occurs, 
+Compression, —Tension, pounds per square inch 


Arch Elevation, 
number feet 


(1) 


Upstream 

Downstream 
Upstream 1,289 (Maximum) 
Downstream —290 
Upstream 1,066 
Downstream 
Upstream 850 
Downstream -47 
Upstream 576 
Downstream —105 


Upstream 
Downstream 
Upstream 
Downstream 229 
Upstream 
Downstream 353 
Upstream —62 
Downstream 307 
Upstream 
Downstream 217 
Upstream 
Downstream 373 (Maximum) 


STRESS (PSI) 


TIME (SECONDS) 


CANTILEVER NO.6 UPSTREAM FACE 
500 
-500 


STRESS (PSI!) 


TIME (SECONDS) 


FIG. 8.—Stress Time History Due Class Earthquake 
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Cantilever Number 
(7) (8) 
(a) Arch Stresses 
287 270 
360 492 213 
481 268 
320 478 398 
426 
280 256 
414 
240 
298 
(6) Cantilever Stresses 
405 
360 110 320 
320 116 378 
—267 
280 
283 
240 
207 
195 
ARCH NO.6 UPSTREAM FACE 
EL. 405 i 
1000 
500 
-500 
2 
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time history the crown near the crest. Typically, the stresses were 
flexural nature with joint openings taking place intermittently indicated 
the tensile stresses. Compressive stresses were well below the ultimate concrete 
compressive strength 4,800 psi (33, 100 kPa), and the region high compression 
localized the crown section the upper half (see Table 7). Joint openings 
took place the upper arches indicated the tensile stresses, shown 
Table The duration these joint openings were approx 0.08 sec. Complete 
separation joints was not indicated. the contraction joints opened and 
closed, arch and cantilever stresses were redistributed varying degrees 


TABLE Earthquake Stresses Time When Maximum Tension Occurs, 
+Compression, —Tension, pounds per square inch 


Upstream (Maximum) 
Downstream 531 
Upstream —283 
Downstream 472 
Upstream 333 
Downstream 268 
Upstream 675 
Downstream 130 
Upstream 552 
Downstream 


Upstream 
Downstream 
Upstream 383 
Downstream (Maximum) 
Upstream 454 
Downstream —298 
Upstream 262 
Downstream 
Upstream 139 
Downstream 157 
Upstream 
Downstream 345 


the arches changed stiffness. The state-of-the-art the analysis arch 
dams unfortunately does not permit exact modeling the nonlinear response 
the dam due the dam’s changing elastic characteristics during the duration 
earthquake. 

order evaluate the effect the reduced stiffness the upper arches 
the redistribution stresses, additional analyses were made adjusting the 
stiffness the central portion the upper arch 25% that used the 
previous analysis. expected, the stresses the upper arches were reduced 
and additional load was taken the cantilevers. The increased cantilever stresses 
were well below the ultimate values allowed. Based these analyses, was 


Elevation, Cantilever Number 
(1) (2) (7) (8) 
(a) Arch Stresses 
405 505 173 
—249 0 
360 531 219 
255 
320 496 
110 480 
280 311 
227 
240 
216 
(6) Cantilever Stresses 
405 
360 172 
—88 160 
320 
382 
280 
210 
240 
126 
195 
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concluded that the dam can withstand the MCE and will function without causing 
uncontrolled releases water from the reservoir. Minor repairs the dam, 
such spalled concrete, should anticipated. 


Summary 


Construction arch dam active seismic zone required careful 
assessment the hazards associated with seismicity. The geology the region 
was reviewed and assessment was made the hazards posed the presence 
three major faults. The hazards were quantified and expressed ground 
motions. The dam was analyzed for two magnitudes earthquakes, called 
Class and Class II. Analyses indicated that the vertical contraction joints 
partially open during vibration, which results redistribution stresses within 
the structure. The results showed that the dam will remain functional and not 
disrupt normal operation after Class earthquake, and after Class earthquake 
minor damage such spalling concrete should expected, but the dam 


will remain intact and uncontrolled release water from the reservoir will 
not occur. 
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INTRODUCTION 


The transverse mixing neutrally buoyant solute that discharged into 
open-channel flow has been the subject many investigations. For uniform 
flow straight, prismatic channels, the advective-diffusion equation describing 
the depth-averaged concentration under steady-state conditions can simplified 
only two terms, expressing the balance between the longitudinal advection 
and the transverse dispersion the solute. For this type flow, the variance 
the concentration-transverse distance distribution increases linearly with 
downstream distance, direct proportion the dispersion coefficient. Thus, 
the dispersion coefficient can evaluated from measured concentration distribu- 
tions. However, when the velocity and depth are not uniform when the 
channel not straight, there transport solute the transverse velocities. 
This transverse advection must included the advective-diffusion equation 
and must accounted for when the dispersion coefficient evaluated from 
experimental data. Taking this into consideration, Holley, al. (5) introduced 
the generalized change-of-moment method for evaluating the dispersion coeffi- 
cient from second moments the distributions solute flux. Their method 
was derived from the advective-diffusion equation written rectangular Cartesian 
coordinates. 

Yotsukura and Sayre (13) generalized this approach writing the continuity 
and advective-diffusion equations general orthogonal curvilinear coordinate 
system which the x-axis follows the meander the river. Then, using 
the cumulative discharge replace the transverse distance coordinate, they 
showed that the equation for steady, two-dimensional mixing can again 
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simplified include only two terms, similar the equation for straight, prismatic 
channels. The advantage this approach that the effects changes 
channel width and depth well channel curvature can readily included 
making particularly suitable for natural streams. These 
effects are much more difficult handle rectangular Cartesian coordinates 
are used. 

This paper presents numerical model based the equation derived 
Yotsukura and Sayre (13). The applicability the model tested using data 
from field experiment. Simulations using different assumptions for the diffusion 
factor are compared. Finally, the dependence the transverse mixing coefficient 
stream variables investigated using published values for natural streams. 


When the solute has mixed fairly well over the flow depth, the two-dimensional 
continuity and advective-diffusion equations orthogonal curvilinear coordi- 
nate system can written 


m 


z 


which time; longitudinal and transverse distance coordinates 
the orthogonal, curvilinear system; depth-averaged velocity components 
the directions; local flow depth; depth-averaged concentration; 
turbulent mixing dispersion coefficients the directions; and 
metric coefficients for the coordinate system. Details for the derivation 
Eqs. and have been given Ref. 13. this study, the x-axis has 
been chosen coincide with the right bank (looking downstream), and 
measured across the stream, orthogonal the x-axis (see Fig. 1). 

When the channel discharge and the solute inflow rate are both steady, the 
time dependent terms can dropped. The longitudinal dispersion term has 
very little influence the transverse mixing under this condition and can also 
neglected (8). Eqs. and are then reduced 


\m, 


Introducing variable called the cumulative discharge which 
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and using place the variable the advective-diffusion equation (Eq. 
transformed into 


the equation presented Yotsukura and Sayre (13) and can shown 
the general form the stream tube models derived earlier papers 

the total discharge the channel, Eq. can written 


LONGITUDINAL 
COORDINATE LINES 


TRANSVERSE COORDINATE LINES 


FIG. 1.—Curvilinear Coordinate System 


well the term which reflects local changes depth, velocity, 
and stream curvature, assumed constant, then Eq. can 
solved analytically. However, for accurate simulations, necessary use 
local values for and solve Eq. numerically. will shown later, different 
approximations for can used, depending the accuracy desired. 


implicit finite-difference scheme proposed Stone and Brian (11) 
solve the advective-diffusion type equations was adopted for the present work. 
When the right hand side Eq. expanded and terms are rearranged, 
takes the form which was analyzed Stone and Brian for numerical oscillations 
and dispersion. The form the equation is: 


called the factor diffusion. includes the turbulent mixing coefficient 
RIGHT 
BANK 
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(i+1,j+1) 


FIG. System for the Numerical Scheme 


yields the following expressions for and 


Referring the grid system shown Fig. and using the discretization 
procedure Stone and Brian, the finite-difference analogue Eq. can 
written follows: 


this scheme, the derivative approximated using weighting coefficients 
€/2, and The other first order derivative approximated using 


and These coefficients have satisfy the following conditions: 


€ 
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weighting coefficients are used for approximating the second derivative 
Crank-Nicolson type approximation used discretize it. 

Stone and Brian have shown that the preceding numerical scheme yields 
solutions closer the analytical solutions (for the case constant and 
when the weighting coefficients take the following optimum values: 


a=b=d= 


They recommended the use the preceding values for the weighting coefficients 
for the general case when and are not constants. 
When the concentration distribution initial station, say known, 
then the concentration distributicn the next station, can calculated 
using Eq. for specific boundary condition the sidewalls. For no-flux 


boundary condition the sidewalls, Eq. can expanded and written 
matrix notation follows: 


in 
S 


(An)? 


a 


When calculating the terms the derivative with respect 
evaluated the same manner the derivative with respect The 
coefficients the tridiagonal matrix can evaluated from the known values 
and and the weighting coefficients. The system equations then 
solved using the Gauss-Seidal method (10) successive substitutions. 

check the accuracy the numerical scheme, the numerical solution 
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was compared with analytical solution which can obtained when 
constant. The initial condition used was concentration distribution 
There were detectable differences between the two solutions. 


Description Test 


field experiment was carried out the Grand River, just south Kitchener, 
Ontario collect data for use with the model. The study reach was 1,500-m 
long with two bends the lower portions (Fig. 3). The streambed was primarily 
sand and gravel with some large rocks. Measurements were made six transects 
which were spaced 250-m apart along the right bank. The stream characteristics 
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FIG. 3.—Test Reach, Grand River Below Kitchener, Ontario 


measured these six transects are given Table The average width, depth, 
and velocity the reach were 59.2 0.506 and 0.353 m/s respectively. 

The tracer used was Rhodamine dye which was injected into small 
tributary wide) located about 230 upstream from the first transect. 
2.5% solution Rhodamine was injected about 100 ml/min using 
constant-rate diaphram pump. The tributary has negligible momentum compared 
the main flow. The dye was mixed throughout the tributary and, entering 
the river, formed small line source the right bank. 

Turner model 111 fluorometer was used the first transect monitor 
the dye concentration. After making sure that steady-state conditions had been 
reached, measurements dye concentration, velocity, and depth were made 
3-m intervals across the transect, using marked cable that was stretched 
across the stream guide. Velocity measurements were taken 0.2 and 
0.8 depths using Ott meter. The fluorometer readings were adjusted for 
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background concentration and also temperature corrected. After completing all 
the measurements the first transect, the field party moved the next transect 


downstream and the whole operation was repeated. Similar measurements were 
made all six transects. 


Using the measurements depth and velocity, the discharge through each 
3-m wide segment the stream was calculated. From these, the cumulative 
discharge with cross-stream distance was calculated each transect. Next, 
the flow distribution curves versus and versus were constructed. 
Using these curves, the variations and with all the transects 
were obtained. The only other piece information required for the solution 

Values for the metric coefficient were obtained graphically. The metric 
coefficient equal the ratio between distances measured along the 
longitudinal coordinate lines and those measured along the x-axis. Therefore, 


TABLE 1.—Hydraulic Data 


Distance Mean Mean shear 


from Channel Average velocity, velocity, U,, 
Transect source, width, depth, meters meters 
number meters meters meters per second per second 
(1) (2) (3) (4) (5) (6) 


definition, 1.0 along the x-axis (right bank). For convenience, the 
longitudinal coordinate lines were chosen coincide with lines constant 
1:1,200 map the river reach, the points constant were marked 
intervals 0.1 along every transect. Curves constant were drawn through 
these points. Distances were then measured along all the lines constant 
and the ratios with distances along the x-axis were computed. this manner, 
the variation with each transect was obtained. 

The dye flux each transect was calculated and compared with the dye 
flux the injection point. For the first five transects, the dye recovery ratio 
varied from For the last transect, the recovery ratio dropped 
77%. apparent reason for this change can found. The dye concentration 
readings were normalized dividing the recovery ratio that the flux 
each transect was equal the input flux. 

Eq. was solved with the numerical scheme described previously, using the 
curve obtained for the first transect downstream the discharge the 
initial condition. Concentration versus distributions for the downstream 


213 49.3 0.47 0.46 0.067 
463 60.0 0.57 0.29 0.074 
713 60.0 0.42 0.42 0.063 
963 55.7 0.28 0.70 0.052 
1213 57.0 0.65 0.29 0.079 
1463 79.5 0.55 0.23 0.072 
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transects were then obtained simulation. These were converted concentration 
versus tranverse distance plots using the curves and were compared with 
the measured concentration distributions. 

The factor diffusion Eq. the product two terms: and 
The term the transverse mixing coefficient which accounts for 
the turbulent diffusive transport well the transport caused differential 
advection. The term accounts for the local variations velocity, depth, 
and curvature. Different assumptions approximations can used for both 
these terms. The term can assumed constant for the whole 
reach can allowed vary from transect transect. One can use local 
values every point use average for the transect 
average value for the whole reach. solving Eq. these approximations were 
tried that their effects the simulated results could compared. The 
different forms which were used are listed the following: 


constant for transect but allowed vary with each transect and 
varies with and 

constant, for the whole reach and varies the variation 


constant for the reach and constant for transect but varies for 
different transects. 


for the whole reach are used for both (uh? 
and e,. The term constant throughout. 

(aU term assumed vary across the stream 
the product the local depth and the mean shear velocity Local values 
for are used. 


Simulations were carried out for all five cases listed the foregoing. 
case which different values for were tried find 
the value which gave the best match between simulated and measured concentra- 
tions the first transect simulated. The solution then proceeded the 
next transect when different values for were again tried. this manner, 
value for was obtained for every transect. The values for varied from 
0.013 The simulated concentrations are compared with the 
measured ones Fig. The agreement can considered very good. 

case the transverse dispersion coefficient was assumed constant, 
for the whole reach. The value chosen for was 0.009 which was 
obtained from averaging the values found case The simulated 
concentrations are also shown Fig. 4(a). can seen that, except for 
the first transect which seems require larger value for e,, the agreement 
between the simulations and the data very good and not too different 
from that for case 

the first two cases, the transverse dispersion coefficieat, e,, remained 
constant across the stream. Case was used investigate the effect allowing 
vary across the stream. this case, was assumed equal 
with being constant for the whole reach. Different values for 
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were tried and the value which gave the closest fit was 0.26. The simulated 
results, shown Fig. are very much the same those for case 
which was assumed constant everywhere throughout the reach. 


MEASURED CONCENTRATION 
PROFILE - 213 m DOWNSTREAM 
FROM SOURCE 
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FIG. 4(a).—Comparison Simulated and Measured Dye Concentration Profiles (Case 
Case and Case 


The last two simulations which were tried involved using average values 
(uh? m,) instead the local values. case average value 
was obtained for each transect plotting the curves versus and 
obtaining the averages graphically. These values were then averaged obtain 
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average value for the reach, for use case The simulated 
profiles for these two cases, for the sake clarity, are plotted Fig. 
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MEASURED CONCENTRATION 
(uh?m,) 

ANALYTICAL SOLUTION 
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FIG. Simulated and Measured Dye Concentration Profiles (Case 
Case and Analytical Solution) 


which separate but identical plot Fig. 4(a). can seen from Fig. 
4(b) that the resulting profiles from these two simulations are very much the 
same. Their agreement with the measured profiles can considered 
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satisfactory, far comparison with field data concerned. However, the 
agreement not good the three cases shown Fig. 4(a), especially 
the last two transects. 


Comparison ANALYTICAL 


the diffusion factor assumed constant for the whole reach, 
Eq. can solved analytically. When the tracer enters the stream from 
tributary, with the present case, appropriate assume the boundary 
solution for this case 


which concentration after the tracer completely mixed across the 
distance. 

Concentration profiles for the various transects were calculated using Eq. 


16. The value for was the same the one used case i.e., 
3.99 ppb and was taken 0.1. These profiles are shown Fig. 
4(b) with the simulations from case and case Overall the agreement 
the analytical solution with the measured data just good those case 
and case actually gives better agreement for the last three transects 
but, the second transect and the first transect, the errors are quite large. 

Basically, the analytical solution and case are the same thing. However, 
case used the measured profile the first transect initial condition 
while this condition could not accommodated the analytical solution. 
This why the two solutions Fig. 4(b) are not identical. 


The results the simulations show that the cumulative-discharge model can 
indeed simulate the transverse mixing natural stream quite satisfactorily. 
The river reach being studied had fairly large transverse velocities some 
sections. This very evident from the fact that between 713 and 
963 the dye plume actually decreased width from about m-21 
Even 713 the plume was still not wide 713 All 
the cases tested, including the analytical solution, were able predict the 
changes the size the plume quite well. one were use the conventional 
advective-diffusion equation such Eq. one would have measure all 
the transverse velocities. these are not taken into account, one can easily 
obtain negative values for the dispersion coefficient. 

From comparison all the solutions, appears that variations the 
dispersion coefficient the downstream direction have bigger influence 
the accuracy the simulations than the possible variations across the stream. 
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The best results were obtained using local values together with 
transverse dispersion coefficient which was constant for transect but 
could vary with each transect. However, allowing vary across the stream 
the local depth did not produce results which were any better than when 
average value was used for the whole reach. 

When average values (transect average reach average) were 
used instead their local values, the simulation results were definitely less 
accurate, especially the transects furthest downstream. Thus, appears that, 
for accurate simulations, one must incorporate the local variations 
and m,. This conclusion different from that Yotsukura and Cobb (12) 
who performed numerical experiments using constant value well varying 
distributions for and found little difference the results. possible 
cause this discrepancy that the distributions used Yotsukura 
and Cobb were symmetrical about 0.5, being maximum the middle and 
zero the edges. The line source was located near the middle and the 
concentration distribution was also nearly symmetrical. Therefore, the lowest 
and highest occurred either where the concentration was very small 
where the concentration gradient was very small and thus their effects 
flux concentration were very small. Where the concentration gradient was 
large, the distribution had value close the average. Under those 
situations, was possible that varying distribution would give much 
the same results constant value for However, the present study 
with the source the bank, the area high concentration gradient occurs 
mainly where the value lower than the average and the results 
are seen significantly different. 

The analytical solution using constant diffusion factor produced results 
which were just good the numerical simulations using average values 
except very close the source. This another advantage 
the cumulative-discharge model—it allows one use analytical solution 
obtain reasonable predictions even for natural streams. 

Thus, for accurate simulations the transverse mixing natural streams, 
one should use local values together with average transverse 
dispersion coefficient which can vary the downstream direction. one 
only interested the size the plume and approximate values the 
concentrations, the analytical solution Eq. may sufficient. 


order use numerical model such the one described this paper 
predict concentration distributions natural stream, necessary 
know the value the transverse dispersion coefficient e,. course, one 
can perform dye test, the value can always obtained from the dye 
data. However, this not always feasible and some method required for 
chosing average for given reach using only bulk stream variables. 

Many measurements transverse dispersion have been made laboratory 
channels. Most the data had been summarized the writers (6) who studied 
the dependence the dimensionless dispersion coefficient the bulk hydraulic 
variables. was shown that the dimensionless dispersion coefficient, was 
function two variables: the friction factor and the aspect ratio W/H. 
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TABLE 2.—Summary Field Data for Transverse Dispersion Coefficient 


Dis- 
persion 
Shear coeffi- 
veloc- cient, 
ity, e,, 
meters squared 
per per 
Data source second second 
(1) (5) (7) 
Yotsukura and 
Cobb (12) Mis- 
souri River near 
Blair 
Yotsukura and 
Cobb (12) South 
River 
Yotsukura and 
Cobb (12) Aristo 
Feeder Canal 
Yotsukura and 
Cobb (12) 
Bernado Con- 
veyance Chan- 
Beltaos (1), Ath- 
abasca below 
Beltaos (1), Atha- 
basca River 
below Ath- 
abasca 
Beltaos (1), North 
Saskatchewan 
River below Ed- 
monton 1.8 x 107° 
Beltaos (1), Bow 
Beltaos (1), 
Beaver River 
Sayre and Yeh (9) 
Missouri River 
below Cooper 
Generation Sta- 
Lau and Krish- 
anppan (6) 
Grand River 


Approximate values, channels are either straight have very small curvatures. 
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Sin- 
uos- 
e,/ ity 
(9) (10) 
0.50 | 1.1 
0.29 | 1.0* 
0.22 | 1.0° 
0.30 | 1.0* 
0.75 | 1.0° 
0.41 | 1.2 
0.25 | 1.0° 
0.61) 1.1 
1.00} 1.3 
3.30} 2.1 
0.26} 1.1 
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the channel width and the average depth. The coefficient can 
take several different forms, such e,/UW, was 
found that the data could best correlated when was plotted against 
W/H. The effect was small and all the data collapsed one curve. How- 
ever, this curve not too useful for selecting values for natural streams 
because most the laboratory data were for W/H less than while for 
natural streams the W/H values are generally much higher. 

was decided investigate the same type dependence could found 
for field data e,. reviewing the published values for natural streams, 
only those which the effects depth-averaged transverse advection, had 
been separated from were considered. quite clear that, for natural streams, 
values which are obtained without accounting for such net transverse advection 


MEAN CURVE 

THROUGH LABORATORY 
CHANNEL DATA, REF 6 


FIG. 5.—Variation with W/H 


are very unreliable. The narrowing the dye plume the present test 
evidence that fact. Therefore, the data which are summarized Table 
are the ones which have been derived from either the generalized change-of- 
moment method (5) some kind stream-tube model similar the one presented 
here. Only data points were found, including the one from the present study. 
The various forms are calculated and are shown Table with other 
pertinent data. 

Fig. plotted against W/H. The mean curve through the 
laboratory channel data summarized the writers also shown. can 
seen that, with the exception two points, the field data are distributed about 
acurve somewhat above the laboratory data. There apparently little dependence 
which varies from 0.015-0.375, and the scatter the data reasonably 
small. However, the two points which not conform with the rest have values 
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which are three times and nine times large those given the mean 
curve. Some explanation this behavior required. 

straight laboratory channels, represents the diffusive transport due 
turbulence well the dispersion resulting from differential advective transport 
due secondary currents which are driven the variations bed shear 
across the channel. natural streams, also has account for the dispersive 
transport due the secondary circulations which are caused the stream 
curvature. natural streams inevitably have some degree curvature, 
reasonable expect that the transverse dispersion due differential advection 
for natural streams will higher, and the values will somewhat 
larger than straight laboratory channels. Thus, quite possible that 
not governed and W/H alone but that additional variable characterizing 
the curvature the particular river reach required. difficult say 
what the best variable use. For wide stream rounding long bend 
which the centrifugal force driving the secondary circulation can assumed 
equilibrium with the pressure forces and turbulent shear stresses, the 
radial velocity component can shown inversely proportional the 
radius curvature the bend. One may thus tempted use the radius 
curvature the additional variable. However, river bends are generally 
not long enough regular enough for the equilibrium conditions established. 
more important factor that river reach would often consist straight 
stretches together with bends different radii curvature. Therefore the 
radius curvature does not seem practical variable characterize 
stretch natural stream. One possible variable which relatively easy 


obtain the sinuosity which can defined the ratio between the thalweg 
length and the down valley distance (7). Using the symbol denote sinuosity, 
the equation for may written 


The values the sinuosity are also listed Table Because lack 
complete information, the values for the Aristo Feeder Canal and the Bernardo 
Conveyance Channel are listed 1.0. However, even though these channels 
have straight sides, the thalwegs are not necessarily straight and for the Aristo 
Feeder Canal has been reported meandering from side side (4). Thus, 
likely that their values sinuosity are slightly larger than 1.0. The same 
holds true for the South River. The values for the Missouri River near 
Blair and the Missouri River near Cooper Nuclear Station (9) are calculated 
using measurements from map. Because the cross sections are not known, 
the channel length taken the thalweg length. Therefore, the actual 
values may slightly higher. 

Inspection the values sinuosity shows that the two points Fig. 
which not conform the mean curve indeed have the largest values 
sinuosity. The one with the largest has the largest sinuosity 2.1. 
The other one has equal 1.3. Ali the other points which are distributed 
about the mean curve have sinuosity ~1.0 1.1. One point has value 
1.2. Although this evidence cannot regarded conclusive, this good 
indication that some variable which characterizes the stream curvature, such 
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the sinuosity, important determining the value natural streams. 
More field data are required, especially for streams with large values sinuosity, 
before any more conclusions can made. There are some published data obtained 
from sinuous laboratory channels but these have very small W/H ratios (mostly 
<10) and would not very useful for natural streams. 

against W/H. There discernable dependence W/H. The data 
can considered distributed about 0.5, although the scatter 
larger than the curve for e,/U,W. Again the same two points with the 
high values sinuosity not conform with the other points. 

There is, present, method for predicting the values for natural 
streams. Therefore, for streams which are relatively free large bends 
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FIG. 6.—Variation with W/H 


other structures such groins which can increase the mixing, the mean curve 
Fig. should give good estimate 


Summary ano 


numerical model the transverse mixing neutrally buoyant solute 
natural channels has been presented. The model equation based the 
stream-tube cumulative discharge concept presented Yotsukura and Sayre 
(13) and can account for changes river geometry and curvature quite easily. 
The validity this model for use natural streams has been verified using 
data from field test. Simulations the dye concentration various transects 
downstream from the source were carried out, using different assumptions for 
the factor diffusion. The simulations all reproduced the narrowing and widening 
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the concentration profiles very well. For accurate simulation the profiles, 
best use local values the velocity, depth, and metric coefficient 
together with transverse dispersion coefficient which can vary the downstream 
direction. However, even analytical solution the model equation can give 


good approximate estimate how the concentration profiles will change along 
the reach. 


Published data the transverse dispersion coefficient for natural streams 
have been summarized. For these data points, the effects depth-averaged 
transverse advection have been accounted for and are not lumped with the 
dispersion coefficient. The dimensionless dispersion coefficient appears 
depend primarily the width-depth ratio W/H and does not seem 
governed the friction factor. All the data points except two seem fall 
reasonably well one curve. suggested that another variable which describes 
the stream curvature, such the sinuosity, also governs the dimensionless 


dispersion coefficient and responsible for the two points not falling near 
the curve. 


For streams which have small curvatures, the curve presented Fig. can 
used obtain reasonable estimate the transverse dispersion coefficient. 
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The following symbols are used this paper: 


a,b 


weighting coefficients; 

depth-averaged concentration; 

fully mixed concentration; 

factor diffusion; 

weighting coefficient; 

dimensionless transverse dispersion coefficient; 
turbulent dispersion coefficient and directions; 
cross-sectional average value e,; 
reach-average value 

friction factor; 

weighting coefficient; 

average flow depth; 

local flow depth; 

weighting coefficient; 

metric coefficients; 

total discharge; 

cumulative discharge; 

sinuosity; 

time; 

mean shear velocity; 

depth-averaged velocities directions; 
width channel; 

width line source; 

longitudinal coordinate; 

transverse coordinate; 

constant; 

weighting coefficient; 

dimensionless cumulative discharge; 
weighting coefficient; 

dimensionless longitudinal distance; and 
function. 
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OVERLAND FLOW FROM TIME-DISTRIBUTED RAINFALL 


INTRODUCTION 


The study hydrology small agricultural and urban watersheds has 
progressed along different lines than that for medium large size watersheds. 
The approach has been more physically based analysis infiltration, overland 
flow, and channel routing. Current techniques for flow analysis are often based 
kinematic wave theory, whereas formerly, they were based the overland 
flow equation derived Horton (6). Horton’s equation might considered 
quasiunsteady analysis because, although used the same basic equations 
those for kinematic wave theory, assumed that the water surface profile 
was always geometrically similar the equilibrium profile. 

Detailed design procedures for small basins have been developed using both 
analytical methods. Hicks (5) and Hathaway (4) based their work the Horton 
equation. Kinematic wave procedures have been described Overton and Tsay 
(12) and Chen and These approaches use rainfall-intensity-frequency 
distributions establish the rainfall input. The average intensity for the 
time-of-concentration used rainfall constant intensity with the Horton 
the kinematic wave flow equations. Using this approach, the rainfall intensity 
decreases time-of-concentration increases; thus, correct time necessary. 
Hathaway (4) showed how rainfall duration influences peak discharge. Time-of- 
concentration depends watershed characteristics, but also varies with rainfall 
intensity which turn varies with time, time-of-concentration determined 
trial and error. Ragan and Duru (13) presented nomograph aid solving 
the kinematic wave relation. Kerby (8) empirically obtained equation based 
results computed Hathaway that apparently appropriate for particular 
rainfall frequency. 


Engr., United States Dept. Agric., Science and Education Administration, 
Agric. Research North Central Region, Watershed Research Unit, 207 Business Loop 
East, Columbia, Mo. 65201. 

Note.—Discussion open until July 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication November 14, 
1979. This paper part the Journal the Hydraulics Division, Proceedings the 
American Society Civil Engineers, Vol. 107, No. HY2, February, 1981. 
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These previous efforts have used rainfall constant intensity estimate 
the overland flow hydrograph. Because constant intensity rainfall unlikely, 
interest determine the influence time distribution rainfall 
the hydrograph overland flow. The kinematic wave equations for turbulent 
flow across plane, impermeable surface are solved using time-varying rainfall 
appropriate for thunderstorms. Peak discharge shown function 
surface length, total precipitation, storm duration, and the time equilibrium 
for equivalent rainfall constant intensity. 


The flow per unit width across impermeable plane surface described 
the continuity equation 


which flow depth (L); discharge per unit width rainfall 
intensity time since start rainfall (7); and distance from 
top plane (L). The total momentum equation difficult apply, 


approximations are usually used. The kinematic wave approach employs the 
relation 


which constants for wide shallow flows. detailed analytical study 


Woolhiser and Ligget (15) led the conclusion that the kinematic wave 
approach appropriate for most formulations. particular, they showed that 
the kinematic wave number greater than 10, the approximation very 
good. The kinematic wave number, given 


which length plane; normal depth outfall with maximum 
discharge; slope plane; and Froude number outfall with maximum 
discharge. This was extended SL/H, Morris and Woolhiser (10). 

The flow the plane may laminar, turbulent, both. reanalysis 
the experimental data obtained Izzard (7) was carried out Morgalli (9) 
confirming this situation. Work Foster, al. (3), indicated, however, that 
use constant Manning’s Chezy’s gave good results when compared 
with their experimental results. the material that follows the flow will 
assumed fully turbulent and that Chezy’s constant. Thus, the constants 
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the solution which may determined using the variable such that, 
the chain rule differentiation, 


— 


and 


To 


Eq. yields the relation the depth flow, rainfall intensity, 
along coordinate The coordinate related the physical time and space 
must established for rainfall intensity. 


THUNDERSTORM 


result extensive studies thunderstorms, the United States Weather 
Bureau (14) published average time distributions for rainfall within high-in- 


1.01" TO 2.00" 


2.01" TO 3.00" <—3.01" TO 4.00" 


4.01" TO 5.00" 


PERCENTAGE OF TOTAL RAINFALL 


PERCENTAGE OF STORM DURATION 


FIG. Distribution Rainfall 1-h Storms [Points Determined Using Eq. 
with 1.37 and 0.37; Curves from United States Weather Bureau 
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tensity storms (Fig. 1). The relation 


storm duration (7); and constants can used describe the curves. 
The value 1.37 and 0.37 fit the 2.01 in./h curve quite 
well, indicated the points Fig. This relationship will used 
describe storms all short duration. 


The depth flow determined integration Eq. 


and are made coincide, the integral the accumulated precipitation 
time, less the accumulated precipitation time, 


Te 


During the storm precipitation given Eq. 


avg 


T 


The location coordinate given Eq. can now determined. The depth 


at 


oO D D 
which can expressed 


HY2 TIME-DISTRIBUTED RAINFALL 
Integration Eq. results 


™%/D 


Evaluating the limits and rearranging terms results 


and describe the water surface profile and hydrograph for the space-time 
coordinates given and during the rainstorm. 

When the rainfall ceases and Eqs. and zero. This indicates 
that the depth, and discharge, longer vary with The space-time 
coordinates this constant flow depth are given 


Examination the results more convenient the equations are nondimen- 
sionalized terms the time equilibrium for equivalent steady rainfall. 
this form, the influence time distribution can compared with constant 
intensity approximation. 

total amount rain, falls time, that the equivalent intensity 


and t>D 
(20) 
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The time equilibrium for this steady rainfall [see, e.g., Eagleson (2) 
Overton and Meadows (11)] 


which the length the plane. The kinematic approximation for the 
momentum equation (Eq. can divided LP/D obtain 


3/2 2/3 pi/3 3/2 3/2 
using the relation for Eq. 21, this becomes 


The quantity sometimes called the runoff intensity, the flow rate divided 
the watershed area. The quantity, P/D, the equilibrium runoff intensity 
for rainfall constant intensity. 

The time coordinate given 


and the space coordinate 


a?/? 3/2 
) 


b+— 


TIME-DISTRIBUTED RAINFALL 


during the rainstorm. After rainfall stops, the space coordinate given 


This form uses the storm duration for the time scale factor, the watershed 
length, for the length scale factor, total fluid length scale 
factor, and the time equilibrium for the equivalent uniform rainfall represent 


the physical characteristics the watershed. 


FIG. 2.—Relation between X/L and for 


The outfall hydrograph determined using the preceding equations the 
following manner. location near the outfall with location selected. 
the start rainfall, the time ¢,, and thus zero. The point indicated 
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determined, using Eq. 25, until the end the plane, reached 
(unless rainfall stops before this achieved). This indicated point 
Fig. The time associated with point used determine the stage using 
the nondimensional form Eq. 13, which 


and the nondimensional discharge, determined using Eq. 23. 
curve C-D. The time which this curve reaches the time 


FIG. 3.—Discharge Hydrographs for Storms Various Durations 


which the total plane contributing the outflow, termed, herein, the 
time-of-concentration. the example depicted Fig. this time occurs before 
the end rainfall. Another situation depicted the curve E-F-G which 
starts with and 0.2. That is, the curve starts after the storm 
has been progress for 20% the duration. The curve point determined 
using Eq. 25. Projecting point accomplished using Eq. 26, however, 
rainfall has stopped for this portion. The stage determined point 
using Eq. and the discharge using Eq. 23. 

Rainstorms that stop before the time-of-concentration achieved, have flat 
peaks. This results because the curve starting +/D will intersect 
the line before reaching Considering the water surface 
profile indicates that the depth flow the right the C-D 
equivalent constant, thus, the flat peak. 


at, 
2.5 
5 
2.0 
1.5 
7 P/O 
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Discharge hydrographs the outfall the plane are shown Fig. Storms 
long duration are depicted the hydrograph for The total watershed 
contributing before the end the storm. The peak occurs the time-of-con- 
centration. The hydrograph then recedes, due the diminution rainfall 
intensity. the end rainfall, D/t, second recession limb appears 
detention storage forms the sole source runoff. 

Storms short duration are indicated, Fig. D/t, 0.5. this 
case, the storm ends before the total watershed contributing the discharge. 
The peak the discharge the end the storm. The hydrograph remains 
that discharge until the water that started leaves the plane, 
after which the falling limb appears. 


10 


FIG. 4.—Time for Total Watershed Contribution Function Time Equilibrium 
for Uniform Rainfall 


The storm duration and the time-of-concentration are equal when 
about 0.87. That is, the time distribution rainfall causes the watershed 
respond more rapidly than does the same amount rainfall distributed uniformly 
over the same time period. The time-of-concentration depends upon the storm 
duration and the physical characteristics the watershed. This time given 
Eq. with x,) and the results which are shown 
Fig. explicit solution cannot determined. The time-of-concentration 
for thunderstorm rainfall always less than for rainfall constant intensity. 
This results from the high rainfall intensity early the thunderstorm coupled 
with the nonlinear flow relation. 

The peak discharge varies with shown Fig. The time distribution 
rainfall has influence the peak discharge the rainfall duration 
less than the time-of-concentration. For the range time, 0.87 D/t, 
1.08, the uniform rainfall leads peak greater than the thunderstorm distribution. 
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This results because the slower time for total watershed contribution with uniform 
rainfall allows more water accumulate. Storms longer duration, however, 
yield higher peaks when distributed according the thunderstorm pattern. This 
THUNDERS TORM jy 
de 1.0 fo 
UNIFORM 
INTENS I TY 
0.1 
0.1 0.5 1.0 5.0 10 
FIG. 5.—Peak Discharge for Storms Various Durations 
because the intensity the early portion the thunderstorm distribution 
greater than that used for the equivalent uniform rainfall. The peak discharge 
determined only this early portion. 
APPLICATION 
The object this effort was evaluate the influence time distribution 
rainfall estimates peak discharge. Although this can inferred from 
Consider watershed with length 500 (152.4 m). The parameter 
in. (5.1 cm) thunderstorm rainfall The time equilibrium for 
this intensity determined using Eq. 21: 
(500)?/* 
(3,600) 
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max P 


The object this study was determine the influence time distribution 
rainfall peak discharge. Usual design procedures use rainfall constant 
intensity for duration equal the time equilibrium. This analysis indicates 
that the peak discharge the design condition will slightly greater for rainfall 
constant intensity than for rainfall with thunderstorm time distribution. 
correct value the time equilibrium essential, however, for estimation 
peak discharge. Difficulty estimating this time parameter limits any value 
considering the time distribution rainfall design applications. 

The constant intensity approximation valid for rainfall durations approxi- 
mately equal the time equilibrium less. the relative duration increases, 
the approximation becomes less valid. Thus, erroneous value for time 
equilibrium can lead significant under design. 
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The following symbols are used this paper: 


rainfall intensity parameter; 
rainfall intensity parameter; 
Chezy constant; 

duration; 


Froude number outfall with maximum discharge; 
depth outfall with maximum discharge; 
intensity; 
parameter kinematic wave equation; 
per unit width; 
peak discharge per unit width; 
radius; 
land slope; 
time; 
flow velocity; 
spatial location; 
variable used integration kinematic wave equation; and 


origin characteristic curve. 
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manuscripts. 
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GRAVEL ANALYSIS FROM 
Closure John 


thank Posey for his discussion and agree that the method need not 
restricted river bed gravels. For the sizing the stockpiles: (1) The 
photograph should taken normal the surface the pile; (2) the surface 
area photographed needs representative sizes within the pile (note 
that rocks the front pile are commonly larger than the average); and 
(3) intermediate axes should measured from the photographs reduce 
the importance rock shape. These cautionary points should also considered 
during other applications the method. 


The author presents mathematical technique for more effective application 
basic hydrologic method commonly used for sizing detention/retention 
facilities. The method sometimes referred the modified rational method 
and various ways applying are discussed Refs. and 

The modified rational method and variations it, such that presented 
the author, are subject potentially serious limitation. The method 
always used with the implicit explicit assumption that the period rainfall 
averaging used the method equal the duration storm. The author 
explicitly makes and states this assumption. the writer’s opinion, this 
not technically correct use intensity-duration-frequency relationships and, 
more importantly, may lead significant underestimates required stormwater 
storage. 


“October, 1979, John Adams (Proc. Paper 14908). 
Fellow, Earth Physics Branch, Observatory Crescent, Ottawa, Canada, 
1980, Kit Burton (Proc. Paper 15229). 


Resources Services, Donohue Associates, Inc., 600 Larry Court, P.O. 
Box 256, Waukesha, Wisc. 53187. 
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Intensity-duration-frequency relationships provide average rainfall intensity 
for periods rainfall and not for discrete storms. For example, assume that 
one enters 10-yr recurrence interval intensity-duration curve for time period 
min and reads intensity 3.0 in./h. This means that once the 
average every (or with 10% probability any given year) there 
will min period continuous rainfall during which the average intensity 
will 3.0 in./h more. does not mean that there will 30-min storm 
with average intensity 3.0 in./h more. Although this 30-min storm 
possible, very unlikely. The most probable situation one which 
the 30-min period during which rainfall averages 3.0 in./h more will 
immediately preceded and followed additional rainfall. 

The assumption that the rainfall averaging time equal the duration 
the storm may have serious implications when rational method concepts are 
used calculate the required volume detention/retention facility. Because 
the method does not account for rainfall that occurs immediately before and 
after the rainfall averaging period, follows that the corresponding volume 
rainfall runoff associated with the preceding and subsequent rainfall not 
reflected the calculated storage obtained with this method. Therefore, this 
method likely underestimate the required volume storage for the specified 
intended design recurrence interval. 

The extent which storage volume likely underestimated not 
generally known. However, the writer’s experience applying the modified 
rational method and comparing the results other sizing methods suggests 
that significant underestimation storage may occur. Therefore, potential 
users the modified rational method, and specific techniques for applying 
the method such that presented the author, are cautioned carefully 
evaluate the results obtained with this technique. 


Baker, R., Water Detention Basin Design for Small Drainage 
Public Works, Mar., 1977. 
Pamphlet—Detailed Steps for Determining Allowable Release Rate and 
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Errata 


The following corrections should made the original paper: 


Page 387, Eq. 18, left hand side: Should read instead 


Page 387, paragraph lines and Should read velocity 


Page 388, Caption Fig. Delete the last equation 1.0 
2 2 


Page 393, Eq. 32, left hand side: Should read instead 
Page 393, Eq. 33a, right hand side: Should read 


l+r 


l+r 


Jump” 


Errata 


The following corrections should made the original paper: 


Page 357, Eq. 11: Should read 1.0 instead 


Page 363, Eq. 22: Should read “‘y, instead 
1980, Tamotsu Takahashi (Proc. Paper 15245). 
March, 1980, John McCorquodale and Abdelkawi Khalifa (Proc. Paper 15250). 
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Discussion Ryszard 


The procedure estimating minimal values the fluctuating pressure causing 
cavitation the stilling basin bottom, presented the author seems over- 
simplified. The author assumes hydrostatic distribution the mean static 
pressure while results the writer’s experiments indicate that pressure values 
the bottom the stilling basin closely the sluice exceed the value 
gy,, (sometimes even 20%). 


FIG. 3.—Pressure Stilling Basin Bottom 


What deserves, however, even more serious objections the author’s very 
estimation the influence fluctuating pressures, upon the cavitational 
process. this writer’s opinion, the duration instantaneous pressure below 
the level should also taken into consideration. According Johnson 
(7), the appearance bubbles the process cavitation requires the time 
question. The mean value the duration instantaneous pressure below the 
level (Fig. may obtained from the following formula: 


“April, 1980, Rangaswami Narayanan (Proc. Paper 15305). 
Dept. Reclamation, Agric. Academy, Poland. 
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Py 


which f(p)dp the intermittency factor (see Fig. 4); and P(p,) 
the mean number runs below the level p,, per unit time. 
For stationary Gaussian functions 


P P 


which the root mean square value dp/dt; the root mean 
square value fluctuating pressure and the mean frequency pressure 
fluctuations. Naturally, the cavitation will appear 


PDF 


FIG. 4.—Skewed Distribution for Pressure 


worth mentioning here that the probability distribution pressure 
Gaussian (normal) and usually positively skewed (As, Ref. 1). This positive 
skewness cannot remain without influence the intermittency factor since 
the case consideration the intermittency should lower than that calculated 
the author for As, Assuming 


easy demonstrate that As, always positive the probability distribution 
the fluctuating velocity (u’) Gaussian: 


Finally, the writer would like draw the reader’s attention calculational 


errors the example provided the author since for 15, m/s 


the values y,, y,, and should respectively 1.63 33.79 6.71 
and 1.84. 


Johnson, E., Journal the Hydraulics Division, ASCE 
Vol. 89, No. HY3, Proc. Paper 3530, May, 1963, pp. 251-275. 
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Associate Members, ASCE 


comparing computed and measured data, one commonly obsessed 
the idea that there some kind discrepancy between the two, the mathematical 
model must adjusted calibrated that finally agrees with the laboratory 
field observations. The method characteristics, particular, has been 
known rigorous, nearly exact scheme for computing open channel flow 
the absence discontinuities, but limited use when applied natural 
streams and channels complex geometry. Its real value lies the fact that 
can serve the touchstone for testing the accuracy and consistency other 
more approximate models, which however, are more tolerant and reliable when 
applied real life problems. For prismatic channel rectangular cross section, 
the characteristic model ideally suited for solving the Saint-Venant equations 
and the writers believe that introduction empirical parameters arbitrary 
assumptions the method not desirable. fact, experience shows that 
best results with respect both measured data and self-consistency the 
model are obtained when the purity the method preserved. 

During the numerous tests the five dambreak models developed Strelkoff 
and the writers Ref. 11, was believed that measurements the Waterways 
Experiment Station the smooth channel were the easiest reproduce, primarily 
because the rigorous estimation the hydraulic resistance. The characteristics 
model, XTICSDAM, ona variable grid using computation nodes with maximum 
node spacing and maximum time step sec, predicted the rate 
advance the surge the dry channel within the WES data 
corresponding the first series. The stage hydrographs also agreed favorably 
with the experimental data for these runs. all cases, the volume error was 
less than 0.5% until recession began, increasing thereafter due 
the neglected volume the trailing edge. 

The author has assumed existing base flow his simulation the dambreak 
floodwave dry channel with apparent success. The generality such fictitious 
initial and boundary conditions questionable, however, the light evidence 
that numerical experiments using progressively smaller base flows initial 
conditions not produce results converging towards the dry bed solution. 
the presence channel resistance, the extremely small depth computed 
the shock front, due downstream depths that are close zero, results 
very large values for the term the governing equations and fact 
renders the entire model meaningless, unless unreasonably small values 

“April, 1980, Cheng-Lung Chen (Proc. Paper 15324). 


Asst. Prof. Water Resources and Civ. Engrg., Univ. Louisville, Louisville, Ky. 
40208. 


Asst. Prof., Dept. Civ. Engrg., Univ. Utah, Salt Lake City, Utah. 
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are used carry out the computations. Whitham’s assumption (12) effect 
removes this singularity quite successfully and, addition, eliminates the need 
for having guess the initial conditions. Apparently, there exists small 
hypothetical value base flow which produces results agreement with the 
dry bed computations and the experimental data depending, however, dam 
height channel geometry and resistance characteristics. Therefore, one would 
have determine empirically the appropriate base flow size for each case, 
which reduces substantially the predictive ability the method. 

addition, the writers would like express their concern with regard 
the use arbitrary base flow and Manning given Eq. for WES 
run 1.2. The resistance characteristics this case are shown Fig. and 
the test was intended under dry bed conditions. However, water trapped 
between the angles used for roughening the surface, created still-water 


MEASURED, WES 1.2 
WET BED, n= 0.05 
DRY BED, n=0.05 
ORY BED, fn from Fig.! 
WET BED, fn from Fig.! 


80 Lele) 120 140 160 180 


DISTANCE FROM DAM, ft. 


FIG. 14.—Advance Trajectories for WES Case 1.2 under Various Channel Conditions 


flow’’ condition average depth 0.06 ft. Advance curves are 
shown Fig. for both dry and wet bed conditions corresponding two 
different modes evaluating channel resistance. Note that the small depths 
encountered near the wave front are way out the range Fig. (or Eq. 
5), which after all was obtained under conditions steady uniform flow. The 
computed advance trajectories indicate that satisfactory agreement with the 
experimental data can obtained one uses existing depth condition 
0.06 downstream the dam and constant Manning value 0.05 (which 
corresponds depth 0.24 that most representative the depths 
encountered over the entire length the wave). 

Another interesting point arises from observation Figs. and the 
negative wave the depth and velocity profiles approach the horizontal tangentially, 
which course incorrect, since the negative wave front represents 
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discontinuity the slope forward characteristics emanating from the undis- 
turbed zone the reservoir. The depression wave front therefore discontinuity 
the slope the depth and velocity profiles too, and interpolation across 
strictly prohibited. The so-called Fig. may very 
well have been the light penalty for this violation. 

the other end the computational domain, the oscillations the velocity 
profiles and hydrographs may simply the consequence the abrupt change 
grid spacing near the wave front associated with the method characteristics 
fixed grid, which results numerical diffraction-like effect the interface 
between computational cells differing too much size. The message 
these that variable, moving grid should used, especially 
under conditions high resistance (11). 


STATION 200 


t= 100 sec 


COMPUTED, STRELKOFF ETA 
© MEASURED, WES 5./ 


280 
DISTANCE, ft. 


FIG. 15.—Comparison Computed and Measured Results for WES Case 5.1 


For the simulation partial breach conditions the author has introduced 
arbitrary parameter, determined trial and error, which besides reducing 
the method characteristics semi-empirical model, unnecessary. Physical 
observation and two-dimensional computations show that the so-called 
veyance has physical meaning the flow through partially breached 
dam (11). fact, its use together with some questionable internal boundary 
conditions, have led the author the conclusion that the characteristic model 
fails the size the breach decreases, which contradicts the most basic 
principles open channel flow. Strelkoff al., (11) have shown that conservation 
mass and energy the breach section are the appropriate conditions for 
the free-flowing breach. solution exists only both sides the breach are 
initially occupied simple waves, i.e., regions uniform flow, and the breach 
itself treated discontinuity depth and velocity. Fig. shows the 
results XTICSDAM for WES case 5.1, corresponding the author’s computa- 
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tions Fig. The hydrograph station 200 (which inference the original 
experimental data appears just upstream the breach), satisfactory 
agreement with the measured data. profile shown the downstream side 
the breach order eliminate the author’s doubts about the accuracy 
the WES data, which, incidentally, remain constant for the period between 
sec and 100 sec. The hydrograph shown Fig. (STA 200) 37% too 
low compared the high side the discontinuity 400% too high compared 
the low side. The fundamental assumption the free flowing breach 
that the flow passes through critical there all times. During the early stages 
dambreak this true for even complete breach, but thereafter, due slope 
and resistance effects, the critical section not necessarily the dam site. 
Imposition critical conditions the breach section for nearly complete failure 
undoubtedly error, since the breach will not small enough choke 
the upstream flow and create the assumed critical conditions. For smaller breach 
widths, however, this assumption will more nearly appropriate for all times 
and thus, the computational results should agree more closely with experimental 
measurements. 

For submerged breach conditions, the author, although results are presented, 
suggests that discharge and depth both sides the 
are the appropriate internal boundary conditions. Across the breach the depth 
may not continuous but should computed judicious application 
conservation energy and momentum (when applicable), coupled with the 
appropriate characteristic equations. addition, one should distinguish between 
conditions submergence the supercritical flow only the downstream 
side the breach discontinuity, which still permits critical section the 
breach, and conditions full submergence, where the upstream flow completely 
coupled with the downstream computations. Results under both these circum- 
stances have been presented Ref. and are further explored paper 
published this journal the writers. 


Discussion Darrel Temple,* ASCE 


The authors are commended for their development straightforward 
approach the complex problem turbulence prediction and analysis. The 
results presented (described the authors preliminary) show this approach 
deserving consideration future investigations. However, their model 
contains inconsistency its assumptions that should recognized and 

1980, Ruh-Ming Li, James Schall, and Daryl Simons (Proc. Paper 
15347). 


Hydr. Engr., United States Dept. Agr., Science and Education Administra- 
tion, Agricultural Research, Southern Region, Water Conservation Structures Lab., Room 
227, Agr. Hall, P.O. Box 551, Stillwater, Okla. 74074. 
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corrected future applications. This inconsistency result using Eq. 
evaluate the viscous term Eq. 10. 

and are statements Prandtl’s velocity defect 
The assumptions classically used analytically derive this velocity distribution 
are presented Schlicting (9) and may summarized as: (1) The viscous 
stress negligible; (2) Prandtl’s mixing-length hypothesis applicable and the 
mixing length proportional the distance from the solid boundary; and (3) 
the shear stress constant over the entire flow depth. 

Since the region applicability Eq. limited that portion the 
flow field where the viscous stress much less than the turbulent stress, its 
use evaluate the viscous term Eq. is, best, questionable. Also, 
the first assumption applied Eq. the result linear variation shear 
stress with rather than the constant required the third assumption. 


TABLE Viscous Term Computed Value 


Variable 


0.001 
0.333 0.026 0.0025 0.0003 


The error resulting from these inconsistencies will depend the relative 
magnitudes and Table constructed using Eq. for “‘typical 
u,, and for flume data, shows the influence the computed viscous 
term small over most the flow depth, suggesting that the results presented 
the authors are valid. evident, however, that the use Prandtl’s logarithmic 
velocity profile prevents the turbulent stress, computed Eq. 10, from satisfying 

Since this behavior primarily consequence applying the first 
results resembling those presented Table would obtained 
for other logarithmic mean velocity profiles derived varying the second and 
third assumptions. 

that portion the flow field near the solid boundary excluded from 
consideration (as must the case when logarithmic mean velocity profile 
assumed), the turbulence-prediction model may made consistent dropping 
the viscous term entirely from Eq. 10. The means computing the time-averaged 
velocity profile then takes secondary importance, although maintenance 
consistent set assumptions desirable. the primary interest 
turbulence-parameter prediction near the solid boundary, the solution not 
clear. Possibilities warranting further investigation include two layer, time- 
averaged velocity distribution models with the viscous term handled consistently 
within each layer, and empirical curve fit models that force the first derivative 
fit the boundary conditions, well forcing the time-averaged 
velocity profile reasonably predicted. 


= 
Data 
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There tendency part hydrology literature consider the conceptual 
flood routing models inherently combined with the finite difference method 
solution. manipulating with the parameters difference schemes one 
can achieve effects absent original models. The author introduced parameter, 
(Courant number) combining conceptual parameter (K) and numerical 
parameter (Ax). Depending the value the parameter, the model 
single linear reservoir produces amplification (for 2), attenuation (for 
2), conservation (for 2). According the author, for the case 
simulate the physical diffusion However, the need for numerical 
introduction diffusion effects the model cascade linear reservoirs 
questionable. The model itself fair approximation the diffusion analogy 
(this statement can illustrated the compatibility impulse responses 
both models). 

The finite difference approach not the only method solution differential 
equation conceptual model. using the impulse response function 
cascade linear reservoirs (Eq. one can conveniently determine the response 
the model arbitrary input signal, making use the following well 
known integral operator 


0 0 

which input signal; output signal; impulse response; and 

term responsible for the impact initial condition upon the output signal 

the time, 

For the case initially relaxed system Eq. reduced the form 
convolution equation [i.e., for all 

the routing example considered the author one can, making use 
integral operators, determine the model responses analytically. The pulse response 
the cascade linear reservoirs can found integrating the impulse 
response the model 


0 


“May, 1980, Victor Miguel Ponce (Proc. Paper 15393). 
Div., Inst. Geophysics, Polish Academy Sci., Warszawa, Poland. 
Hydro., Inst. Geophysics, Polish Academy Sci., Warszawa, Poland. 
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which duration time rectangular pulse input signal; and unit 
step response the model (S-hydrograph). The values the integrals from 
Eq. can determined analytically for integer values the parameter 
with the help the following equation (12). 

x"e“dx =—x"e" x" dx 
The theoretical solutions prove that the model acts diffusionlike manner, 
whereas amplification conservation the amplitude are not inherent properties 
the model and result from the numerical method finite differences. 

the general case real value, and arbitrary input signal, analytic 
evaluation the integral from Eq. not possible. However even then 
simulation can easily performed with the help pocket-sized calculator. 
Study more sophisticated convolution models solved with the help 
programmable calculator was presented Keefer (15). 

Integral operators enable calculation the model response for any number 
linear reservoirs directly from the input signal. One does not have calculate 
the transformation the flow intermediary reservoirs. 

Assuming the kernel function properly defined, routing convolution 
always stable. Therefore the choice the time interval, not constrained 
any stability condition. The time step, At, warranting the proper accuracy 
the numerical convolution routine depends upon the dynamics the input 
signal and the system impulse response function. 

the convolution integral evaluated the following discrete representation 


then the value the time step (or the analogate the Courant number 
At/K) should sufficiently small ensure that all discrete moments 
approximate well the moments respective continuous time functions. 
particular, the continuity condition 


i=l 


i=l 
and the condition compatibility the time step chosen and the dynamics 
the input signal 


i=l 


= > | x(t’) dt’ 
i=l 0 


should fulfilled. 


easy see that large values are not concern convolution 
modelling, since such case the bulk the impulse response function (or 
range which variations input signal occur) may skipped over. 

The formula for the parameter (reach length corresponding linear 
reservoir, i.e., the so-called characteristic reach length) was obtained Ponce 
the Muskingum model. The physical sense the parameter can found 
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more directly by: (1) Simplifying the hydrodynamical equation flow (14); 
and (2) comparison impulse responses the conceptual model cascade 
linear reservoir and the linearized Saint Venant equations [moment matching 
method 

Another point that the writers’ opinion deserves discussion the author’s 
attitude with respect off-centering the difference schemes for the Muskingum 
model and for the kinematic wave approach. the Muskingum model has 
approximate the kinematic wave model one has use the value parameter 
equal 0.5 (17). 

will shown later cascade Muskingum models with parameter 
0.5 gives better approximation the kinematic wave model. Let apportion 
the overall modeled reach the length into subreaches, each the length 
AL. Let consider the cascade Muskingum models, each having the 
parameters 0.5 and AL/c L/(nc), which the constant 
wave celerity. tends infinity, the cascade Muskingum models tends 
the kinematic wave model. This can clearly demonstrated means 


the transfer function approach. The transfer function for multiple Muskingum 
model reads 


1-— 


can proved with the help the theorem, this transfer 
function tends for the transfer function the kinematic wave model 


shown Ref. the values the parameter different from 0.5 are 
the consequence hydrodynamics; thus, the off-centering the scheme for 
the Muskingum model physically interpretable. 
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Dooge, I., Theory Hydrologic Technical Bulletin No. 
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The authors discuss the theory the linear hydrograph method and correctly 
state its limitations. They point out that many urban watersheds are within 
these limitations. The writer has found that there are many practical applications 
for similar method (17) published 1975 the analysis and design urban 
and some rural drainage systems. 

The comparative advantage single linear reservoir model (SLR) used 
the Santa Barbara Urban Hydrograph (SBUH) method the ease application 
and the few parameters required. the SBUH method, calculated 
the time concentration the watershed, the and 
therefore most civil engineers are already familiar with the procedure. The time 
concentration most practical applications can estimated from the physical 
characteristics the watershed and, unlike 15, independent 
storm characteristics. The writer believes that the assumptions necessary 
apply Eq. real watersheds, which not have planar surfaces constant 
rainfall intensities, are more accurate than approximations used estimating 
times concentration. 

The SLR model and SBUH method require little computational effort compared 
with other hydrograph methods. Routing rainfall excess through linear reservoir 
does not require computer, and can even done graphically (20). step 
step procedure for calculating hydrograph the SBUH method follows: 


The storm rainfall amount and distribution are determined. The storm 
divided into time periods each with duration the time interval 
selected. The appropriate amount rainfall assigned each time period 

The impervious portion the watershed (Imp) which effectively hydrauli- 
cally connected estimated. 

The pervious infiltration rate inches per time increment (/) determined. 

The watershed area acres (A) and the time concentration hours 
(Tc) are determined. 

Runoff depths for each time period are calculated: 


Pervious area runoff, >=0 
Impervious area runoff, (26) 


“May, 1980, John Pedersen, John Peters, and Otto Helweg (Proc. Paper 
15430). 


“Engr.-Mgr., Santa Barbara County Flood Control and Water Conservation District, 
123 Anapamu St., Santa Barbara, Calif. 93101. 
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Total runoff depth for each time period, R(t) 


The Hydrograph”’ ordinates cubic feet per second 
for each time period are obtained multiplying the runoff depths for each 
time period the drainage area and conversion factor: 


1.008 


I(t) 


The final hydrograph ordinates cubic feet per second are obtained 
routing the instantaneous hydrograph through single linear reservoir with 
time delay equal the time concentration 


Note that Eq. applied only once for each time period. 


inwhich 


addition the author’s and writer’s work, the SLR model used 
the SBUH method has been tested against observed data Florida Golding 
(19) with good results. The writer believes the method appropriate many 
practical cases where complete hydrographs are required, such design 
stormwater retarding basins. Therefore, important emphasize how easy 
for nonspecialized engineer apply that the model may used 


more widely. 


19. Golding, L., Synthesis the HNV-SBUH presented 
the June 19-20, 1980, SWMM Users Group Meeting, held Toronto, Canada. 

20. Wilson, W., Graphical Flood Routing Transactions the American 
Geophysics Vol. 21, Part 1941, pp. 893-898. 


The following symbols are used this discussion: 


drainage area, acres; 

infiltration rate pervious area, inches per time period; 
instantaneous hydrograph ordinate, cubic feet per second; 
hydraulically connected impervious portion watershed decimal; 
rainfall time period, inches; 

final hydrograph ordinate for time period cubic feet per second; 
pervious area runoff, inches; 

impervious area runoff, inches; 

runoff depth for time period, inches; and 

time concentration, hours. 
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Discussion Richard Hawkins” 


The author has innocently credited the originality Eq. Smith (8). 
The basic authorship might more appropriately attributed Cheng-lung Chen, 
who first derived the expression 1975 (9,10). This mixup due large 
part oversights the writer review stages Hjelmfelt’s manuscript. 
While exposed Chen’s 1975 work several times, the writer was not fully 
conscious the presence Eq. and thus overlooked correcting the 
earlier dated citation. fact, several other papers the writer (2,11,12,13) 
between 1975 and the present contain the same equation—or near equivalents— 
without reference Chen’s work. Thus, spirit honest embarrassment, 
apologies are extended here Chen for slighting his earlier contributions. 
Hjelmfelt’s paper has served vehicle promote this recognition. 

Future researchers probing curve number algebra would well study 
Chen’s pioneering work detail. Lacking even earlier dated claim, perhaps 


would appropriate hereafter refer Eq. 4—and closely allied forms—as 
i.e. 


dt’ 


should noted that the foregoing equation correction what 


probably typesetting error the published paper. The exponent should 
refer the entire denominator, and not just 
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DISCUSSION 


The authors have presented several useful models for predicting the steady-state 
performance cooling ponds with different spatial structure. their last 
paragraph they suggest that the general use these models will for the 
estimation average, such monthly, performance conditions. This use will 
provide basis for comparing initial pond designs, design modifications, etc. 

They also suggest that steady-state models can employed with shorter-term 
data approximate the transient response ponds. However, because the 
response pond transient meteorological input governed its depth, 
the appropriate time interval for averaging input data the order 
rather than the residence time implied. fact, assuming linear heat transfer 
and constant station operation, the correct transient solution may viewed 
(6) time series steady-state solutions evaluated with exponentially 
filtered time series T,, using filter parameter approximately 
correct transient response obtained one averages meteorological input data 
over time interval, pc,H/K. The value typically between 
day and week which, coincidentally, the same range typical pond 
residence times given 

also noted (6), that the maximum intake temperature predicted with 
steady-state model can sensitive the averaging interval the 
spatial structure assumed for the pond (e.g., longitudinally dispersive versus 
laterally recirculating, plug flow versus well-mixed, etc.). This suggests that 


the transient response pond important consideration, even the 
design stage. 


Adams, E., and Koussis, A., Analysis for Shallow Cooling 


Journal the Energy Division, ASCE, Vol. 106, Proc. Paper 15737, Oct., 1980, pp. 
141-153. 


1980, Gerhard Jirka and Masataka Watanabe (Proc. Paper 15449) 

and Principle Research Engr., Ralph Parsons Lab. for Water Resources 
and Hydrodynamics, Dept. Civ. Engrg., Bldg. 48-325, Massachusetts Inst. Tech., 
Cambridge, Mass. 02139. 
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